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ABSTRACT 

Star- forming galaxies display a continuous distribution of specific star formation rates (sSFR) which, 
in the simplest approximation, can be reproduced by the superposition of two log-normal distribu- 
tions. The first of these encompasses the main sequence of star-forming galaxies, the second one a 
rarer population of starbursting systems. We show that the sSFR-distribution of starbursts can be 
regarded as the result of a physical process (plausibly major mergers) that takes the mathematical 
form of a log-normal boosting kernel and that enhances star formation activity. In this article we 
explore the utility of splitting the star-forming population into main-sequence and starburst galax- 
ies - an approach we term '2-Star Formation Mode' (2-SFM) framework - for understanding their 
molecular gas properties. We show that variations of star formation efficiency (SFE) and gas fractions 
among star-forming galaxies take a very simple, redshift-independent form, once these quantities are 
normalized to the characteristic value of an average, massive main-sequence galaxy. The change in 
SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected 
for merging events. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane 
that separates more clearly into loci for starbursts and normal galaxies than observed in a plot of 
(s)SFR vs. M*. Starbursts with the largest deviations (>10-fold) from the main sequence, like many 
local ULIRGs, are not average starbursts, but even rarer events with progenitors having larger gas 
fractions than typical main-sequence galaxies. We statistically infer that the gas fractions of typical 
starbursts are reduced by a factor of two to three with respect to their direct main-sequence progeni- 
tors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained 
more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 con- 
version factor aco in the SFR-M* plane and provide strong evidence that the higher sSFR of distant 
galaxies is a direct consequence of the higher gas fraction of these systems. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: spiral - galaxies: ISM - 
surveys 



1. INTRODUCTION 

Studies over the last decade of star-forming galaxies 
(SFGs) have revealed a positive and tight correlation 
between their current star formation rate (SFR) and 
stellar mass M*, which is intimately linked to the 
integral of the preceding star forma tion (SF) activity. 
Initia l ly observed at low redshift (e.g. 
[2001 ISalim et all 120071 IWvder et al 



, Brinchmann et all 
20071 ). this 'star- 



forming main s equence' was soon shown to be present 
out to z ~ 2 (INoeske et all 120071 lElbaz et all 120071 
iDaddi et al.l I2007LJ ) . Subsequent work on the relation 
between SFR and in SFGs charted its evolution 
to z ~ 2.5 usi ng different S FR t r acers and se l ection 
criteria (e.g. iDamen et al.1 [2001 iDunne et al.l I2009L 
Pannella et al.l I2009L iSantini et ah 12009 ~|Kaiisawa et al 



2010L IQliver et al.l I2010L 
20 111 iRodighiero et al.l 



lElbaz et al 



2011 



.. iLee et al 
T20TT[ iWuvts et all 1201 ll 



iWhitaker et al.l I2012D . has explored factors affect- 
ing t he exact shape and dispersion of the m ain sequence 
(e.g. iKarim et all 120111 iSalmi et al.l 120121: T. Bschorr 
et al., in prep.) and has traced it out to even higher 
redsh ifts z ~ 3-4 (e.g. IDaddi et all 120091 iMagdis et all 
|2J310j). The existence of a scaling relation between SFR 
and M* throughout much of cosmic time implies that 
the (mass-dependent) assembly history of SFGs is char- 
acter ized by a high degr e e of homogeneity and simplicity 
(e.g INoeske et all 120071 iBouche et all 120101 iPeng et all 
120101 lLeitnerll2012L iBehroozi et al.1 l2012h and that the 
stron g decline of the co s mic s t ar formation rate d ensity 
(e.g. iRedd-y 



y fc Steidell l2009l IRodighiero et all 12010, 
Karim et all 120111 . iMagnelli et all 120111 iCucciati et all 
2012t and references therein) since z ~ 2 reflects the 

of the majority of the SFG 



* E-mail: mark.sargentOcea.fr 

1 CEA Saclay, DSM/Irfu/Service d'Astrophysique, Orme des 
Merisiers, F-91191 Gif-sur-Yvette Cedex, France 

2 Department of Physics, University of Oxford, Keble Road, 
Oxford OX1 3RH, UK 

3 Smithsonian Astrophysical Observatory, 60 Garden Street, 
Cambridge, MA 02138, USA 

4 Max-Planck-Institut fiir Astronomie, Konigstuhl 17, D- 
69117 Heidelberg, Germany 



and 

uniform SFR-cvolution 
population rather than a decreased frequency of episodic 
starburst events. Nevertheless, the study of starbursts 
remains central to understanding the nature of inter- 
acting galaxies and the physics of merging events which 
may produce the m ost luminous sources at all redshifts 
(|Sargent et al.H2012[) . 

Star formation activity at a rate which locally oc- 
curs only in strong starbursts is common in massive 
main-sequence galaxies in the distance universe. Hence 
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alternative strategies to pure luminosity-selection arc 
required for obtaining a census of bursty SF activ- 
ity at high redshift, e.g. bas ed on their posit i on in 
the SFR vs M± p lane (e.g. iRodighiero et al.1 1201 lL 
Whitaker et al.ll2012D or based on the i r mor phology (e.g. 
Kartaltepe et al.ll2012l iKaviraj et all 120131 ). The latter 



20081 [2010a| ITacconi et all l20l"ol IGeach et all I20TH 
Tacconi et al.l I2013D . If, as discussed in iRenaud et all 



approach relies on the assumption that starbursts are 
generally triggered by interactions between galaxies, as 
observed in the lo cal universe (e.g. iSanders et afl 119881 
iBarton et aHteOOOt ). while main-sequence galaxies would 
represent a 'normal', secular channel of stellar mass 
growth in galaxies that is f ueled by the steady accretion 
of cold, primordial gas (e.g. lBouche et al.1l2010D . 
A series of studies on the star-forming population have 
improved our understanding of normal (MS) galaxies 
and starbursts and are in qualitative agreement with 
this picture. Sta rbursting sources a re more compact 
on a verage (e.g. lElbaz et all 1201 lL iRuiopakarn et al.l 
[20Tll) than main-sequence galaxies which have a stellar 
stru cture that is well described by exp onential disks 
(e.g. lWuvts eTaI|[201lT ISalmi et al.ll2012f ). They display 
deficits in the intensity of infrared (IR) spec tral features 
(e.g. PAH bands or in the far-IR [Cllj-line; lElbaz et all 
120 111 iGracia-Carpio et alJl20rl . and have warmer IR 
spect r al energy distributions (SEDs; e.g.lHeisler fc Vaderl 
1991 ISanders fe Mirabel U9M IChapman et al.1 12001 



Elbaz et al.l 120111 iMagdis et al.l I2011L iBethermin et al l 
20121 and references therein). These are telltale features 



of intense and spatially concentrated star formation as 
is expected to occur in interacting or merging galaxies 
whe re gravitational torques funnel gas to th eir centers 
(e.g. iMihos fe Hernquist]|1996l iHopkins et aXll2006t) . 
The efficiency with which gas is converted into stars in 
such settings may be up to an order of magnitude higher 
than in the extended gas reservoirs which fuel star 
formation activity i n normal galaxies ou t to z ~ 2 (e.g. 
Daddi et al.l I2010U iGenzel et all 120101 ITacconi et all 



This strong contrast in star formation effi- 
ciency (SFE = SFR/Mg as ) is often taken as one of 
the most clear-cut manifestations for the existence of 
two distinct star formation laws - a secular mode in 
main-sequence galaxies and a starburst mode charac- 
terized by short depletion tim escales (<100Myr, e.g. 
iSolomon fc Vanden Bouti[2005t and references therein). 
Whether or not such a bimodalit y represents the physica l 
reality has been questioned (e.g. iNarayanan et al.ll2012| ) 
on the grounds of discrete 'concordance' values being 
assumed for the CO-to-H2 conversion factor aco and, 
secondly, due to the expectation that star formation 
laws at a basic level should be expressed in terms of 
volumetric quantities rather than observationally more 
easily accessible surfa ce densities of SFR and gas (e.g. 
IKrumholz et al.l I2012D . The lack of known sources 
with SFEs between those measured for normal disk 
galaxies and strong starbursts could also be a selection 
effect: initial CO follow-up observations of high-redshift 
galaxies targeted only highly luminous sources experi- 
encing 'bursty' star formation (sub - millimeter galaxies 
and QSOs; e.g. IQmont et al.l Il996l iFrayer et all 119981 
Walter et all 120031 IGreve et al.l I2005L iMaiolino et al.l 



20071 ITacconi et aLll2008l ) and followi ng improvements in 



the sensitivity of millimeter receivers (IChenu et al.l[2007L 
iPerlev et al.l 120 111 ), dedicated studies of t ypical main- 
sequence galaxies were undertaken (e.g. iDaddi et al.l 



(|2012fl . the gas density distribution function - which 
reflects the turbulence-driven structure of the interstel- 
lar medium (ISM) - is a crucial factor in determining 
the shape of star formation laws, then intermediate 
SFEs should indeed occur in, e.g., minor mergers or 
in certain stages of galaxy interactions when the gas 
density distribution is not modified from the steady 
state as strongly as during final coalescence. However, a 
dichotomy in the distribution of SFEs could still occur 
if the timescales for such variations wer e short (e.g. 
iTevssier et aT][2010llBournaud et al.ll2011al) . 
In this paper we consider a large sample of local and 
high- redshift SFGs which we use to extend the '2 Star- 
Formation M o de fr amework' (2-SFM) introduced in 
I Sargent et al.l (|2012t henceforth abbreviated as 'S12') 
to the molecular gas component of SFGs. The 2-SFM 
framework relies on basic observables (e.g. the evolu- 
tion of sSFR in main-sequence galaxies or their stellar 
mass distribution) and correlations between observables 
(e.g. the star-forming main sequence or the Schmidt- 
Kennicutt relation). Our goal is to describe how the star 
formation efficiency (SFE) and the molecular gas content 
of galaxies are related to their location with respect to 
the main sequence, i.e. to their specific star formation 
rate (sSFR), which is the main diagnostic of 'starbursti- 
ness' within the 2-SFM framework. We will postulate 
that there is a different law linking the SFE and sSFR for 
main-sequence galaxies and starbursts, and derive these 
laws from comparison with the observations. We will 
show that the population of massive SFGs that reside on 
the main sequence has similar molecular gas properties 
across a broad range of redshifts (z < 3) and we will use 
our detailed description of the starburst population and 
its SFR-'boosting' developed in Section 14.21 to demon- 
strate how in the 2-SFM framework a bimodal behavior 
in terms of SFE arises naturally even in the absence of 
discrete star formation laws for normal galaxies and star- 
bursts. The description of SFE in the SFG-population 
developed in this paper forms the basis for the predic- 
tion of molecular gas mass functions and CO luminosity 
functions in a companion paper (Sargent et al. in prep.; 
Paper II, henceforth). 

The outline of this manuscript is as follows. Section [2] 
introduces the observational dataset we use and how it 
was homogenized. We then employ this reference sam- 
ple in Section [3] to calibrate galaxy-scale star formation 
laws - both in terms of observables (Lir & L' C q) or in- 
trinsic quantities (SFR & Mh 2 ) - for galaxies at low and 
high redshift. These calibrations depend on the adopted 
CO-to-H2 conversion factor acd and the correspond- 
ing systematics will also be assessed in Section |3] The 
mathematical description of the starburst population in 
the 2-SFM framework, and its relation to main-sequence 
galaxies is the focus of Sect. 0] where we derive the dis- 
tribution of the burst amplitudes - the 'boost-function' 
- that transforms a theoretical population of pure main- 
sequence star forming galaxies into the observed distri- 
bution of sSFR. We discuss what physical mechanisms 
could produce this boost-function and consider in par- 
ticular the possible link between starbursts and galaxy 
mergers. In Section [5] we will combine the redshift- 
independent Schmidt-Kcnnicutt law derived in Section [3] 
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with the evolution of the sSFR distribution from S12 to 
construct prescriptions for the relative variation of molec- 
ular gas properties of normal and starburst galaxies that 
are particularly simple (and self-similar) once they are 
referred to the properties of the average main-sequence 
galaxy. Our results are presented in three main blocks: 
the SFE and gas fractions of main-sequence galaxies and 
starbursts are the subject of Section [5~T1 and [5721 respec- 
tively; Section 15.31 focuses on the CO-to-H 2 conversion 
factor and its variation within the main sequence and 
among starbursting systems. We then discuss and sum- 
marize our findings in Sections |6] and [7] 
Throughout this article we adopt the WMAP-7 
cosmology (f) m = 0.273, + £l m = 1 and H = 70.4 
km s~ 1 Mpc~ 1 ; lLarson et al.ll201lT ) . Star formation rates 
and stellar masses are given for a IChabrierl (I2003T) ini- 
tial mass function^ (IMF). All literature values have 
been adapted accordingl y. Metallicities are given on the 
IKewlev fc Dopital (|200a henceforth 'KD02') scakQ and 
(molecular) gas mass estimates include a 36% correction 
for helium. 

2. DATA 

We discuss two different kind s of data sets in this sec- 
tion. To begin with (Sect. I2.1j) . we describe individual, 
CO-detected SFGs which we will utilize to establish the 
basic scaling relations (e.g. for SFE or gas fractions) that 
link the molecular gas content of massive (M* > 10 10 M©) 
SFGs to fundamental galaxy properties like (s)SFR or 
M*. In Section 12.21 we introduce statistical samples of 
SFGs at z~l and 2. These will subsequently be used 
(1) to visualize/simulate complete samples of galaxies 
that obey the aforementioned scaling relations, and (2) 
to extend the analysis to fainter galaxies where the va- 
lidity of such scaling relations can be verified with image 
stacking. 

2.1. The reference sample of individual star-forming 

galaxies 

Our 'reference sample' of normal galaxies at redshifts 
z < 3 comprises 91 sources from the recent literature (see 
Sects. 12.1.11 & 12.1. 2jt . We complement these with lo- 
cal and high-redshift starbursts with measured CO-to- 
H2 conversion factors (see Sect. I2.1.3|) . These are es- 
sential for a further investigation of the notion that the 
'bimodality' of star for mation is particular l y pronounced 
in te rms of SFE (e.g. iDaddi et al.l l2010at iGenzel et al.l 
l201Gh . 

2.1.1. Normal galaxies: low redshift CO-detections 
The HERACLES survey (|Lerov et al.l [200l [200l 



120131 ) targeted the CO(J=2— >1) transition in nearb 
(£><15Mpc) THINGS galaxies ([Walter et al.l [200^ 
with the IRAM 30 m single-dish telescope. Here we 
select 20 galaxies with spiral galaxy morphology and 
stellar mass M* > 10 10 M Q from the H ERACLE S sam - 
ple. Stellar masses (converted to the IChabrierl ()2003l ) 

Logarithm ic mas ses and SFRs based o n a Salpctcr l|1955l ). 



scale) and morphol ogical information are taken from 
the compilat ions of iSkibba et al.l (|2011l ) or iLerov et al.l 
(|200l [2009h . or from the NASA/IPAC Extragalactic 
Database^ if not listed in either of the former. Metal- 
licity estimates for m ost of the select e d HE RACLES 
spirals are provided in iMoustakas et all (|2010f ). The IR 
(8-1000 //m) luminosities attributed to the HERACLES 
galaxies a r e bas ed on the photometry reported in 
iDale et all (120071) and have been calculated following 
eg. 22 in lDraine' fell (12001 . 

We augment the local main-sequence galaxies from the 
HERACLES data set with a subset of g alaxies from the 
first r elease of the COLD CASS survey (jSaintonge et al.l 
120111 ) for which an accurate IR luminosity could be 
calculated thanks to the presence of a c ounterpart in 
eithe r the IRAS Faint Source Catalog (v2; iMoshir et al.l 
19921) or the AKA fll/Far-Infrared Surveyor (FIS; 



a lKroup a (2001) and a Baldry & Glazcbrook (20Q3) IMF are con- 
verted to the Chabrier scale by adding -0.24 dex, dex and 0.02 dex, 
respectively. 

7 When necessary, metallicity information from the literature 
was con verted to the KD02 cali bration by means of the prescrip- 
tions in IKewle v & EllSor] l[2008l ). 



Kawada et al.l 120071) all-sky survey Bright Source 
Catalog (vl.O: lYamamura et al.ll2010D. An add i tiona l 
cut in stellar mass (taken from iSaintonge et al.1 120 111 ) 
at M* = 1O 1O M excluded all less massive sources. 
Among the 222 sources in the first COLD CASS data 
release 32 fulfill these criteria; they lie in the redshift 
range 0.025 < z < 0.05 and are all late-type galaxies 
with CO(J=1^0) fluxes measured by the IRAM 30 m 
telescope at S/N > 4. IR luminosities for these galaxies 
were computed u sing the spectra l ener gy distribution 
(SED) library of iCharv fc Elbazl pOOll ) and allowing 
renormalization of the templates when fitting the 
reliabk0 IRAS or AKARI photome try at A rfiS t > 30 am i 
(for details on the IR SED-fitting see lHwang et al.ll2010D . 

2.1.2. Normal galaxies: high redshift CO-detections 

CO-transitions in main-sequence galaxies at z 3> have 
been targe ted bvlGeach et all (| 2009l [201lt CO(J=1^0) 
at z ~ 4).lDaddi et al l (l2010al lbl: COf J=2^1) at z -0.5 
& 1.5). iTacconi et al.l (l20ldl~C Of J=3^2) at z-1.2 & 
2.3), and lMagdis eTaLl (|2012al: CO(J=3->2) at z~3) - 
all using the IRAM Plateau de Bure (PdBI) interferom- 
eter - leading to line flux measurements at S/N>4 to- 
ward 32 of 39 observed galaxies (for the remaining galax- 
ies 3 a upper flux limits are available.) The stellar masses 
of these sources are in the range 5 x 10 10 < M*/ M Q < 5 x 
10 11 , as determined by SED fitting of the near-UV to 
near-IR broad-band photometry. 

The derivation of SF Rs varies am o ng th e differ ent afore- 
mentioned studies. iGeach et al.l (|2009l . 120 111 ) estimate 
the IR-luminosity from the flux of the 7.7 /xm emission, 
as co nstrained by the S pitzer IR spectrograph (IRS), 
while [Daddi et all (|2010bD based the luminosity measure- 
ments for their z^ 0.5 sources on Spitzer/MIPS 24 /im 
fluxes. S FR-estimates for zy l.b sBzK galaxies pre- 
sented in Daddi et al.1 (|2010a[ ) are an average of dust 
corrected UV luminosities, mid-IR continuum luminosi- 
ties from 24 /im imaging and VLA 1.4 GHz radio con- 
tinuum fluxes, all of which were found to give consis- 
tent SFR-estimates. iMagdis et al.l (|2012al ) adopted a 
similar averaging approach but were able to add Her- 
schel/PACS and SPIRE photometry to constrain the 

8 http://ned.ipac.caltech.edu 

9 Flux quality flags are either 'high' or 'moderate' for IRAS 
sources and 'high' for AKARI sources. 
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dust-emission of their z ~ 3 LBGs. We adopted all 
these SFR (ira) measurements without further modi- 
fications but in the interest of maximal sample homo- 
geneity ch ose to re-compute L tr values for the z~1.2 
sources in iTacconi et al.l (|2010t) as these authors report 
SFRs that are based on a combination of extinction- 
corrected optical emission lines, UV flux and mid-IR 
photometry. To derive new Ltr estimates we searched 
the Far Infrared Deep Extragalactic Lega cy (FIDEL, PI: 
M. D ickinson, see also data description in iMagnelli et al.l 
[2001 coverage of the Extended Groth Strip for associ- 
ated Spitzer/MIPS 24 fim detections and converted these 
to IR luminosities with the ma i n-sequ ence star-forming 
galaxy IR SED of lElbaz et al.l (120 111). Tw o of the ten 
sources with 1.2 from the ITacconi etall (12010ft sam- 
ple h ad VLA 1.4 GHz d etections from the AEGIS 20 cat- 
alog (jlvison et al. 2007), such that the 24 /im- and radio- 
based SFR co uld be average d . For sources at z^2.3 in 
the sample of ITacconi et al.l (|2010l) . neither IR nor ra- 
dio flux measurements we re available; w e henc e used the 
original SFRs reported bv lTacconi et al.l (|2010l ) which are 
based on extinction-corrected Ha luminosities and which 
we con verted to an IR-luminosity following IKennicutt] 
(|1998bn . 

To summarize, we have compiled a sample of 91 mas- 
sive (M* > 10 10 M Q ), main-sequence galaxies with CO- 
detections of which 57% are local {z < 0.05) systems and 
the remaining 43% redshifted to 0.35 < z < 3.2. 

2.1.3. Starbursting galaxies with measured aco 

Starbursting galaxies with the same information as 
available for our reference sample of n ormal galax- 
ies are listed in the recent study of iMagdis et al.l 
(12012b!; starburst-lik e galaxies at 2.3 < z < 4) and in 
ISolomon etaLl l|1997t local starbursting IRAS ULIRGs). 
An accurate assessment of the behavior of SFE during 
starburst episodes - one of the main aims of this paper - 
requires an observational determination of the CO-to-H2 
conversion factor aco- In the case of the local ULIRGs, 
we thus restrict ourselves to nine objects - VII Zw 31, 
Arp 193, Arp 220, Mrk 273, 00057+4021, 02483+4302, 
10565+2448, 17208-0014, 23365+3604 - with two inde- 
pendent measureme nts of aco: one based on dynami- 
cal constraint^ by iDownes k Solomon! (|1998l) and the 
other based on large velocity gradient (LVG) r adiative 
transfer modeling bvlPa padopoulos et al.l (|2012f ). 
The flux of the CO( J=l— >0) transition toward our sub- 
sample of nine z < 0.07 /iM.S'-detected ULIRGs was 
measured with the IR AM 30 m single-dish telescope by 
ISolomon et al.l (|1997l ). We computed their IR lumi- 
nosities (and thence SFRs) using spectroscopic redshifts 
from NED and all available IRAS photometry, follow- 
ing th e standard recipes provided by Sanders k Mirabell 
p996t their Table 1). Stellar masses - which are particu- 
larly important for the characterization of these systems 
in terms of sSFR, our prime indicator of "starbursti- 
ness" - have been p ublished for some of the ULIRGs 
in our sample (e.g., lElbaz et al.l 120071 Ida Cunha et al.l 

10 In the following we adopt a mass-to-light ratio aco that is 
given by the ratio between gas mass and dynamical mass within a 
~l-3kpc region encompassing both the inner, high-density nuclear 
disk/ring and an outer, lower-density disk (wi th volume filling fac- 
tor ~0 .1 for the gas) of the ULIRGs modeled in lDownes fc Solomon! 
(fT998l ; see their Tables 3 & 9). 



l20T0llHowell et al.l2010llU et alj2012T) . but to our knowl- 
edge no single study has done this consistently for all 
sources of interest. We hence re-est imated stellar masses 
for the IDownes k Solomon! (I1998D ULIRGs based on 
2MASS if-band fluxes (jSkrutskie et al.1 12006! ) and pre- 
scr iptions for ma s s-to-li ght ratios, M+/L k, as derived 
by lArnouts et al.l (|2007t their eq. 2) and Uuneau et al.l 
(|201 It their eq. B 2). The stellar masses we adopt in 
the following for the IDownes k Solomon! (|1998l) ULIRGs 
were obtained by averaging the estimates calculated ac- 
cording to these two prescriptions. They agree well 
with the available literature measurements (median off- 
set 0.1 dex). 

In addition to the nine low-redshift starbursts just de- 
scribed we also include the three high-z sub-millimeter 
galaxies GN20 (z = 4.05), SMMJ2135-0102 (z = 2.325) 
and HERMES J105751. 1+573027 (z = 2.957) in our anal- 
ysis. The re cent determination o f their conversion 
factor aco in IMagdis et al.l (|2012bf ) by means of the 
M gas /Mdu S t-ratio t echnique (see also iLerov et~aU 1201 lL 
IMagdis et ail 12011 relied on: (i) a far-IR SED/dust 
emission that is accurately constrained by Herschel and 
millim etric continuum observations (see IMagdis et al.l 
12012b! for a detailed listing), and (ii) CO(J=1->0) line 
fluxes from the ( J) VLA (for GN20; see lCarilli et"aT1l2010l 
iHodge et al.l 120121) and the Green Bank Telescope (for 
SMMJ2135-0102 and H ERMES J105751.1+5 73027; sec 
ISwinbank et~aT|l2010l and lRiechers et al.ll2011l resp.). All 
three sources have sSFR-cnhancements with respect to 
the main sequence of at least a factor three, as con- 
st rained by the optical to near-IR and IR SED-fitting 
of IMagdis et al.l (|2012bf ). Furthermore their CO-to-H 2 
conversion factors are systematically lower than that of 
the Milky Way, similar to the values typically measured 
in interacting local ULIRGs. 

2.2. Statistical samples of star-forming galaxies in 
GOODS-South 

With the purpose of demonstrating the applicability 
of our recipes for computing molecular gas properties for 
observed galaxy samples, based on individual measure- 
ments of stellar masses and SFRs, we use two samples 
of A'-selected galaxies in the GOODS-S fi e ld, at z~l 
and 2 taken from th e work of iDaddi et al.l (l2007al) and 
iDaddi et al] (|2007bl: see also |Sahmet_aL| [2012| foT more 
details on the z ~ 1 sample) . The same samples were 
used in the rece nt pap ers by IMagdis et al.l (|2012bD and 
IMullanev et~aTI (|2012l) . We refer to the original papers 
for details of how stellar masses were derived, based on 
empirical recipes using colors and absolute luminosities. 
The SFRs of the galaxies at z ~ 1 and 2 are based on 
24 /im and UV observations, respectively, and are known 
to compare well on average with ot her tracers including 
Herschel-based SFR-measurements (jDaddi et al.l I2007al 
lElbaz et al.ll2010llReddv et al.ll2012!) . 

3. GALAXY-SCALE STAR FORMATION LAWS: 
CORRELATIONS BETWEEN SFR AND GAS 
MASS AND THE ASSOCIATED OBSERVABLES 

Re cent reports (e.g. IDaddi et al.ll2010bl iGenzel et al.l 
2010 : of a systematic offset between the power law re- 
lation linking the overall surface density of SFR and 
gas mass (the Schmidt-Kcnnicutt (S-K) law Esfr oc Sg as ; 
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Fig. 1. — : Observed correlation between measures of star formation rate and molecular gas content for massive 
(M* > 10 10 M Q ) main-sequence galaxies and starbursts at low and high redshift. Low redshift sources (filled/open 
black symbols - normal galaxies; grey crosses - starbursting (U)LIRGs) arc from the HERACLES, GOLD GASS and 
IRAS surveys. Rcdshifted galaxies from the literature (see legend) are plotted in color (orange - z ~ 0.4; red - z ~ 0.55 
& z~1.5; blue - z~1.2 & z^2.3; magenta - starbursts at 2.3<z<4), with arrows indicating 3cr upper limits for 
CO non-detections. Open blue an d cyan symbol s indic ate modified SFR estimat es (see section [2 1.2 1 for details) based 
on radio or IR data for sources in iGenzel et al.l (|201QO and iTacconi et al.1 (|2010t) . The shape of the symbol indicates 
which CO transition was detected toward the individual sources (cf. legend in lower right corner of panel b) . 
(a): Correlation between infrared luminosity (£ir) and CO-luminosity (£co(j=i-s.o)' standard excitation corrections 
- e.g. iDannerbauer et all (|2009f ) . ILerov et ahl (|2009f ) - were applied to J > 1 transitions) with the best- fitting relation 
derived for the main-sequence galaxy sample plotted as a solid black line. The strong starbursts considered here (cross 
symbols) are on average offset to higher SFRs by a factor of three (dashed line). 

(6): Inverse, integrated Schmidt-Kcnnicutt relation between SFR and molecular gas mass (A/ mo i ), the latter having 
been derived based on cither (i) observational determinations of aco (available for sources with boxed symbols) or (ii) 
using a metallicity-depcndcnt conversion factor (see text of Sect. I3.2|) . The dispersion about the best-fit linear trend 
(solid black line) for main-sequence galaxies is approximately Gaussian with a dispersion cr(M mo i.) ~ 0.2 dex (red curve 
in inset). Dashed line - offset locus with approx. 15-fold higher SFE for strong starburst galaxies. 



ISchmidtJll95i lK enmeu ttl ll998ah of normal galaxies and 
starbursts were highly influential in shaping the notion 
of 'bimodal' star formation. These findings are subject 
to two systematic uncertainties. Firstly, the measured 
offset between normal and starburst galaxies depends on 
(potentially population-dependent) recipes for CO-to-H2 
conversion factors which are hard to estimate for a sta- 
tistically significant number of SFGs, especially at 0. 
Secondly, the sampling of the S-K plane obtained as a 
result of targeted CO follow-up observations toward se- 
lected SFGs is patchy. Constructing a reliable and sta- 
tistically representative sampling of the distribution of 
SFGs in the S-K plane is thus not only important to ex- 
plore different modes of star formation. The observed 
S-K law is also often referred to as a benchmark for the 
performance /validity of recipes for ISM processes in sim- 
ulatio ns (e.g. lRobertson &: Kravtsovll2008L [Monaco et al.1 
l2012t ) and hence used to gauge our understanding of the 
underlying physics itself. 



In this section we return to our literature compilation of 
low- and high-redshift main-sequence galaxies with CO- 
detections that we presented in Sects. 12.1.11 and 12.1.21 
We re-measure the slope and dispersion of the galaxy- 
scale star formation law and in doing so for the first time 
are able to in corporate ano-measu rements from the re- 
cent study of [Magdis et al. (2012b]) for a fraction of our 
reference sample of normal galaxies. Rather than using 
SFR and gas mass surface densities, we consider the sim- 
pler relations between integrated quantities, namely the 
total SFR and molecular gas mass or the corresponding 



observables, Lir and L' ( 



CO(J=1->0)- 



For the rest of this 



article we will use '£c ' as a shorthand for the line lu- 



minosity L' ( 
12 CO. 



CO(,/=1->0) 



of the first rotational transition of 



3.1. L' 



L 



IR 



V CO(J=1^0) 

We begin with the relation between the observables, 
Lir and L' co , that are the starting point for cstimat- 
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ing SFRs and molecular gas content of SFGs. As de- 
scribed in Sect. 12.11 all galaxies considered in the fol- 
lowing have stellar masses M*>1O 1O M0. Current ob- 
servations of BzK- s elected main-sequenc e gala xies (e.g. 
iDannerbauer et ahl 120091 . lAravena et al.l l2010T l suggest 
that typical excitation corrections for the first two higher 
order transitions J=2 — »1 and 3— » 2 are_r 2 i = 0.8-0.9 and 
r 3 i ~ 0.5. Similarly, iLerov et"aH ()2009|) find an aver- 
age J=2^1/J=1^0 line ratio of 0.8 for HERACLES 
galaxies. In Fig. [TJi we plot the accordingly corrected 
CO(J=l— >0) luminosities of local (black and white sym- 
bols) and redshifted sources (colored symbols) against 
their IR luminosity. We then fitted the CO line luminos- 
ity as a function of IR luminosity. While performing a 
regression of L m on L'qq would be more natural (as rep- 
resenting the relation between cause and effect, i.e. M mo i, 
and SFR, rcsp.) our choice is motivated by the aim to 
provide recipes for the molecular gas content and associ- 
ated tracer emission beginning with the currently obser- 
vationally more e a sily a ccessible SFR-measurements. A 
iBucklev fc JanTesI (fl~979l; ' BJ' hereafter ) regr ession (im- 
plemented as described in llsobe et al.1 119861 ) , which al- 
lows for a statistically correct treatment of the 3 a up- 
per de tecti on limits for six ga laxies from iTacconi et al.l 
(pOlOh and lGeach et~aTl poll , gives 



log 



-^CO(J=l->-0) 

K km/s pc 2 



= a\ 



with (1) 



[ax; Pi) = (0.18 ± 0.02; 0.84±0.03) for normal galaxies. 

The dispersion about this best-fit trend line in the y- 
direction is 0.21 dex. In performing the linear regression 
we have down-weighted sources detected in CO( J=3— >2) 
by a factor of two due to the large excitation corrections 

bl- 
under the assumption that starbursts follow a correla- 
tion with identical slope, we use our reference sample of 
starburst galaxies (see Sect. 12 . 1 .3[) to solve for the nor- 
malization of eq. [T] which best reproduces their average 
offset. We find 



(«i;/3i) = ( 



-0.28^; 0.84) 



i.e. an offset of 0.46 dex or approx. a factor 2.9 with 
respect to the locus of main-sequence galaxies. This 
similar systematic d ifference was already indicated by 
iSolomon et al.l (|1997t see their Fig. 3) in their pioneering 
analysis of CO-emission in nearby ULIRGs. Local main- 
sequence galaxies with the IR luminosities of starbursting 
ULIRGs are, however, exceedingly rare (e.g. S12) such 
that this difference could also have been explained by a 
double power-law nature of the star formation law or a 
single, steeper relation (Ltr oc -L' co 13 ) owing to different 
probability g as density distributions in mergers and nor - 
mal galaxies ((Narayanan et a l. 2008, Jun eau et al.ll2009l ). 
The advent of CO line flux measurements for high-z nor- 
mal galaxies with ULIRG-luminosities has since added 
another piece of evidence in support of a systematic off- 
set (e.g. IGenzel etld1l2ljl0l ) 

3.2. M mol . vs. SFR 

The integrated S-K law linking the molecular gas mass 
(M mo i.) and SFR is expected to have a different slope 
or curvature than the correlation between logarithmic 



luminosities L' co and Ltr unless the average CO-to-Bi2 
conversion factor 



«CO(J=1^0) 



mol. 



CO(J=1->0) 



is a constant. Evidence to the contrary has been 
presented in numerous stu d ies, the most recen t of 
which are ILerov et all (|2Qlll ). iSchruba etall 1)20121 ) and 
IGenzel et al.l (|2012D who show - albeit with large un- 
certainties - how aco decreases with metallicity among 
galaxies in the Local Group and in massive main- 
sequence galaxie s at z < 2.5, respectively. Note that 
ISandstrom et al.l (|2013| ) find a similar trend for decreas- 
ing aco with increasing metallicity when considering 
spatially distinct regions within nearby galaxies. In 
the next paragraphs we discuss a scheme for assigning 
metallicity-dcpendcnt aco _v ahies to the normal SFGs in 
our reference sample (Sect. 13.2. ip and then proceed to 
fit the resulting relation between SFR and M mo i. in Sect. 
13.2.21 We close this section with an assessment of the ro- 
bustness of the galaxy-scale star formation law obtained 
in this way (Sect. I3.2.3p . 

3.2.1. Statistically inferred CO-conversion factors 

For nine of the 91 main-sequence galaxies in our ref- 
erence sample the CO-to-H2 conversion factor has been 
measured and found to be broadly consistent with the 
(negative) power-law relation between aco and metal- 
licity found for galax y samples in the nearby Universe 
(jMagdis et al.ll2012b| ). For the remainder of the popu- 
lation we assume that their aco too follows a relation 
with metallicity Z that has the form 



log(aco) = v + Clog { z /z @ ) 



(2) 



and that metallicity can be deduced in a statisti- 
cal sense from stellar ma ss and SFR as p r opose d by 
IMannucci et al.l (|2010t ) and lLara-Lopez et al.l (|2010|) who 
describe the interdependence of these three quantities as 
a 2D manifold in M*-SFR-Z space (see Fig. Hi). In the 
following we use this 'fun damental metall i city r elation' 
(FMR) as parametrized by IMannucci et aTl (|2010f ). Since 
starbursts constitute a small fraction of the star-forming 
population, the FMR primarily reflects the dependence 
of metallicity on SFR and M* for main-sequence galax- 
ies. As the stellar mass and SFR of the CO-detections 
in our reference sample are known (within observational 
errors) we can use the FMR to statistically infer mctallic- 
ities Z(SFR, M*), and thence CO-conversion factors aco 
for each of the galaxies in our reference sample^. This 
enables us to calibrate the M mo i. vs. SFR relation in- 
dependently from the correlation between L' co and Lir 
derived in Sect. 13.11 Note that more general expectations 
for aco-varia-tions in the M* vs. SFR plane that account 



11 IMannucci ct al. (2010) originally were only able to study 
the FMR at 2>0 for massive (M*/Mq > 10 10 ) field galaxies. 
New work has sin ce extended the relation to lower stellar masses 
llUresci ct al. 201i|) and separately verified it s validity in the cluster 
environment at z~ 1.4 IIMag riiii ct al. 2012). The FMR is gener- 
ally assumed to hold over the range < z < 2.5, i.e. should apply to 
almost all galaxies in our reference sample. Beyond z ~ 2.5 conflict- 
ing evidence for constancy (e.g. Dcssaugcs-Zavadsk v et aTJ 1 20111 . 
IRichard et ahirSOlll, fLara- Lopez et al.H2012j. IBelli et alj|2013) and 
evolution (e.g. ILaskar*^t**aT . 2011, Sommariva ct al. 2012) of the 



FMR has been presented. 
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Fig. 2.—: Overview of recipes used to assign CO-to-H 2 
conversion factors, qco, to observed and modeled galax- 
ies based on their SFR and stellar mass. 

(a) : Funda mental metallicit y rela tion (FMR) as de- 
termined by IMannucci et all (|2010t thick line section) 
and continuation of the FMR assumed in our analytic- 
empirical modelling to regions of parameter space where 
the mathematic a l desc ription of the plane proposed in 
IMannucci et al.l (|2010t) diverges due to the absence of 
data. Note that the extension beyond /Uo.32 € [9.25, 11.4] 
affects only a minority of the galaxies modeled in the 
present work. 

(b) : Recently proposed relations between metallicity 
(expressed in multiples of the solar m etallicity Zq in 
the system of iKewlev fc Dopitaj |2002f) and aco based 
on simulations (iFeldmann et all I2012D o r observations 
(jSchruba et al.l 120121 IGenzel et all 120121 : for the latter 
study the fit derived for both low- and high-z SFGs (light 
grey, dashed) and that for a sample restricted to z > 1 
SFGs (light grey, solid) is shown). All recipes predict a 
Milky Way-like aco at approx. solar metallicity but di- 
verge significantly at lower/higher enrichment due to the 
different measured slopes (see annotations beside trend 
lines) . 



for (1) changes of the conversion factor within both the 
normal and starbursting galaxy population, and (2) the 
relative importance of these to classes of SFGs depend- 
ing on the location in the plane, are the topic of 15.3.21 
Existing literature consistently reports the normaliza- 
tion v of the aco vs. metallicity relation in eq. 
[2] to be such that Milky Way-like conversion factors 
a co = 4.4M Q /(Kkm/spc 2 ) are reached around solar 



metallicity. Measurements and expectations for the slope 
£ span a larger range which we illustrate in Fig. [2Jd. In 
view of these poor constraints we took a different ap- 
proach to estimating slope and normalization of the rela- 
tion between aco and metallicity. The details of this cal- 
culation (and the underlying assumptions) are described 
in Paper II where we use the 2-SFM framework to infer 
the redshift-evolution of the molecular gas mass function 
and CO luminosity function (LF) starting from the stel- 
lar mass function of SFGs. The 2-SFM prediction for the 
z ~ CO LF can be directly compa red to the local CO 
LF measured bv iKeres et al.l (|2003l ). In particular, the 
Schechter function parameters of the observed local CO 
LF are sensitive to the choice of v and £ in eq. [2] The 
best-fit slope and normalization we obtain in this way are 
£~-0.9 and v~ 0.7 (or aco — 5 at solar metallicity Zq) 
and we subsequently adopt these values to assign conver- 
sion factors to all galaxies in our reference sample that 
lack an observational determination of aco- It should 
be clear that our strategy for constraining the shape of 
the aco vs. metallicity relation is subject to certain 
systematics (e.g., the assumed scatter of the correlation 
between aco and metallicity) which we discuss in Paper 
II and which are the reason why we do not provide error 
estimates on v and £. However, our main goal here is to 
determine the integrated S-K relation for massive main- 
sequence galaxies and we will show in Sect 13.2.31 that it 
is quite insensitive to the slope and normalization of eq. 
[2] due to the almost fla t mass-metallicity rel ation in this 
stellar mass range (e.g. lTremonti et al.ll2004f ). Moreover, 
it is interesting that with slope £~-0.9 and normaliza- 
tion v ~ 0.7 the scatter of the S-K relation is smaller than 
with any other of the recipes for aco shown in Fig. [2J). 

3.2.2. Integrated Schmidt-Kennicutt laws 

Multiplication of the excitation-corrected CO- 
luminosities of Fig. [TJt with (i) observed CO-to-H2 
conversion factors if available or (ii) the statistical 
CO-conversion factors discussed in the previous section 
provides a measure of the molecular gas mass for each 
of the main-sequence galaxies in our reference sample. 
These measurements are plotted against their SFR 
in Fig. []J). At the stellar masses considered here, 
the IR-excess -Lir/Luv of mai n-sequence galaxies is 
in general ~10 or larger (e.g. iWhitaker et al.l 120121 : 
M. Pannella et al., in prep.), leading to a nearly 1:1 
correspondence between Lir and SFR, as indicated 
by the lower and upper scale for the x-axis of Fig. 
[H>. We conve rt IR luminosities to SFRs following the 
prescription of iKennicutti (|1998bf ). 

BJ-regression, applied to the data in Fig. [T]d, returns 
a very similar logarithmic slope as for the correlation 
between L' co and Lir due to the relatively narrow range 
of metallicities spanned by our low- and high-z data: 



log 



Mr, 



Mr. 



Ot2, SFR + 02 log 



SFR 

A'%/yr 



with (3) 



(a2,SFR;/32) 



(9.18±0.02; 0.83±0.03) for normal galaxies 



.00 



+0.29. 
-0.10' 



0.83) for strong starbursts. 



Here the line parameters for the starbursts in our refer- 
ence sample were derived by solving for the normalization 
under the assumption of an identical, slightly sub-linear 
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Fig. 3. — : Robustness assessment of the integrated inverse S-K relation calibrated in Fig. [T](see also eq. [3]). (a) Changes 
in the best-fit S-K re lation, depending on the metallicity-dep endence, of the conversion factor, ago, reported by 
recen t observational ([Genzel et al.ll2012l iSchruba et al.ll2012[ ) and numerical ((Narayanan et al.ll2011l iFeldmann et al.1 
120121 ) work. (Filled/hatched areas encompass the 95% confidence region o f the corresponding bes t-fit relations.) 
The additional dependence on the CO surface brightness, Wco, proposed by INaravanan et al.l (|2011[ ) introduces an 
implicit dependence on galaxy size (long dashes - CO-flux averaged over optical half-light radius; short dashes - CO- 
flux averaged over two optical half-light radii) . 

(6) Extension of high-z S-K relations to the highest SFRs (^20M Q /yr) reached by the z^O main-sequence galaxies 
in our calibration sample. (Grey symbols in backgro und reproduc e data p lotted in Fig. [TJb). After inclusion of the 
stacked samples of z ~ 2 main-sequence galaxies from|Magdis e t al.l (j2012bt colored symbols), the z ~ 2 S-K relation is 
sampled over two orders of magnitude in SFR. The stacked samples bridge the gap between local and z > calibration 
sources and are aligned with the best-fit S-K relation determined using these, thereby providing evidence of a universal 
S-K relation for massive main-sequence galaxies. 



slope of the star-formation law for main-sequence and 
starbursting galaxies. Alternatively, in units of £tr, 



log 



M mo i. 



= a 2 , ir + log 



-Hi 



with 



(4) 



f (0.90±0.02; 0.83±0.03) for normal galaxies 
(a 2 , ir; #2 J - | (_ .28+0;29; 0.83) for strong starbursts. 

The dispersion of the correlation is 0.20 dex and almost 
equal for the low- and high-z subsample s as shown in the 
inset p anel of Fig. [Tp where we plot the lKaplan fc Meierl 
(|1958[ ) estimator for the cumulative distribution func- 
tions of the offsets A(M mo i.) of the individual measure- 
ments from the best-fitting trend line in eq. EH The 
average offset of the reference starbursts with respect to 
the locus for normal galaxies is 1.18+ ; 2 gdex or roughly 
a factor f 5 in terms of SFEP^I. It should be emphasized, 
that this bimodality is arbitrary and merely reflects the 
properties of our s mall an d incomplete selection of strong 
starbursts. In Sect 15. 1^21 we propose an empirical descrip- 
tion of the starburst population that allows for a contin- 

12 In this section we ha ve used the dynamical constraints from 
IDownes fc Solomon! ( 119981) to convert L' n n to M mo i for the local 
starburst ULIRGs. As discussed in Sect. 15.3.11 the average offset 
from the integrated S-K law doe s not change for this sampl e when 
we use the aco- vamcs given by Papadopoulos ct al. (2012). 



uous enhancement of SFE, depending on the importance 
of burst-induced star formation activity. 

3.2.3. Robustness of calibrated M mo l. vs. SFR relation 

The integrated S-K law for main-sequence galaxies 
which we calibrated in the previous section (eq. [3} con- 
stitutes a key ingredient for the description of the gaseous 
component of SFGs in the 2-SFM framework. It is thus 
essential to ascertain that the shape of our best-fit S-K 
law is not strongly dependent on assumptions made dur- 
ing the calculation of, e.g., gas masses. 
A first potential cause of systematic uncertainty is our 
choice of the exponent of the relation between aco 
and metallicity discussed in Sect. 13.2.11 In calibrat- 
ing the integrated S-K law, we adopted a dependence 
aco °f Z~ 09 which is intermediate bet ween the ob- 
served trends that range from Z° (e.g. iBolatto et al.l 
12008( 1 to Z~ 2 - 7 (e.g. Ilsraellfl997t l. As representative ex- 
amples of these different measurement^^ we show in Fig. 

13 Observational studies generally parametrize «co _var i a tions 
as a function of an absolute value of oxygen abundance rather than 
relative to the solar metallicity. The inferred oxygen abundances 
may vary significantly dependin g on the metallicity calib ration and 
abundance diagnostic used fe.g. IKewlev fc E llison 2003). In addi- 
tion to these systematic uncertainties concerning the normaliza- 
tion of the metallicity scale, transfo rmations between metal licity 
systems following, e.g., the recipes in lKewlev fc E llison (20(]3) a l so 
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[3^ the star formation laws that we obtain when apply- 
ing the simulation-based recipe of lFeldmann et all (|2012t 
aco oc Z~ ' 5 ), and those observationally determined by 
IGenzel et ail (|2012t a co oc Z" 127 ) and ISchruba et al.l 
aco « 2 2 ), to the normal galaxies in our refer- 
ence sample (where we again have assigned mctallicities 
using the FMR). We find that the slope of the differ- 
ent S-K laws are very similar and that their normaliza- 
tion only varies by 0.2 dex, such that they are all con- 
sistent within uncertainties among each other and also 
with our best-fit star formatio n law as given by eq. [3[ 
We also tested the prescription iNaravanan et al.l f|2011h 
developed based on their simulations of disks and merg- 
ers. These authors parametrize aco as a function of 
metallicity and CO surface brightness, Wco, which in- 
troduces an implicit dependence on galaxy size. Using 
optical size measurements from the literaturj 1 "^! and the 
CO-fluxes of our reference galaxies we calculated CO-to- 
H2 conve rsion factors and t hence molecular gas masses 
following INaravanan et al.l (|2011[) . Our subsequent fit 
to the data showed that - for reasonable assumptions 
about the relative spatial d istribution of optical an d CO 
emission, see Fig. [3^ - the INaravanan et al.l (|2011f) pre- 
scription also leads to a star-formation law that is fully 
consistent with our preferred integrated S-K relation (eq. 
[3]). The good general agreement between all these differ- 
ent recipes is due to the flatness of the mass-mctallicity 
relation at M* > 10 10 M . 

A second source of systematic uncertainty is our as- 
sumption that massive galaxies at all redshifts align 
along a single S-K law. By combining high-redshift 
gala xies from the P HIBB S survey with COLD CASS 
data iTacconi et al.l (j2013l ) recently presented an alter- 
native scenario of parallel and linear S-K laws which 
are c haracterized by an S FE that increases with red- 
shift. ITacconi et al.l (120131 ) based their gas surface mass 
densities on CO( J=3— ^2) fluxes for z > 1 SFGs and uni- 
formly applied a Milky Way-like conversion factor to all 
sources in their sample, regardless of mass, SFR and red- 
shift. Both the strong excitation corrections applied to 
the high-z SFGs which dominate the high-E gas regime of 
the S-K relation and the universal aco are in principle 
uncertain enough to bring about seemingly systematic 
shifts between the high- and low-redshift galaxy popula- 
tion in the S-K plane. Given the fact that there is little 
overlap in Ltr between our own subsamplcs of low- and 



change the curvature of the aco vs - metallicity relation. As an 
example, a log-linear relation b etween ago and metallicity de- 
rived using the R23 indicator (e.g. Pilyugin 2001, Kcwl cv fc Dopital 
120021 ) may become a convex or concave function of metallicity when 
converted directly to a system based on the N2 diagnostic (e.g. 
IDenicolo et~aTi r2002. Pcttini & Pagel 200|). In view of these com- 
plications we chose to renormalize all literature determinations of 
the aco vs - metallicity relation considered here to solar metallicity 
(see Fig. [2p). In terms of oxygen abundance, solar enrichment gen- 
erally corresponds to a value of log(0/H)+12 £3 8.7, but this may 
change somewhat depending o n the metallicity calibra tion (e.g. 
log(Z Q )+12 = 8.9 in the scale of lKewlev fc Dopitall2002T) 

14 Optical half-light radii for main-sequence galaxies in 
our reference sample we r e der i ved u sing the following litera- 
ture sources: ILerov et alj f2 008. 200iJ) f or HERACLES galax ies; 
IForster Schreiber et a" . (2009, 2011]) and IGenzel et all pOlOf l for 
SINS galaxies; IDaddi et al.l H2010al ) for CO-detected sBzK galax- 
ies. No size information was available for the z ~ 0.4-0.6 galaxies 
from Gcach et al. (2011) and Daddi et al. (2010b), nor for galaxies 
in the COLD GASS sample (Saintoiigc ct al. 20n|). 



high-redshift galaxies plotted in Fig. [Tb, we cannot rule 
out a series of offset and conceivably also curved star 
formation laws. By includi ng average SFE-cons traints 
from the stacking analysis of lMagdis et al.l (j2012bl ) , how- 
ever, it is possible to bridge the luminosity-gap between 
z > 2 and local SFGs, as shown in Fig. [3b . In combina- 
tion w ith the BM/BX-selected galaxies of ITacconi et al.l 
(120101) , the z ~ 2 S-K relation thus spans one and a 
half orders of magnitude and is seen to extend contin- 
uously into the parameter space of intermediate-redshift 
(0.4 < z < 0.6) and local disks without evidence for a dis- 
continuity. When considering individual de tections, the 
z~0.4 galaxies of lGeach etahl (120091 120111) are aligned 
wit h the z~ 1.2 and 1.5 sample of ITacconi et al.l (|2010l ) 
and IDaddi et a l. (2010a). Moreover, the skew and scat- 
ter around our universal S-K law are very similar for the 
subsets of z ~ and z 3> galaxies in our reference sam- 
ple (see inset of Fig. [T}d). Based on these observations, 
we conclude that the assumption of a single, slightly sub- 
linear relation between SFR and M mo \, is presently a 
valid working hypothesis. 

4. THE INNER WORKINGS OF THE 2-SFM 
FRAMEWORK: MAIN-SEQUENCE GALAXIES 
AND BOOSTED, STAR-BURSTING SOURCES 

While the S-K relation we found in Sect. [3] for main- 
sequence galaxies is very well defined, it is much less 
obvious which concrete form of the star formation law 
should be used to describe the sparse and scattered star- 
burst data. Our approach to interpreting the incomplete 
information on these sources will be to statistically link 
them to a synthetic and complete starburst population 
where we are able to relate the starburst properties to 
the pre-starburst (main-sequence) state. We describe in 
the following the '2 Star-Formation Mode framework' (2- 
SFM) which is the basis for establishing this link. 

4.1. Basic ingredients and successes of the 2-SFM 
framework 

2-SFM is a simple and self-consistent scheme for the 
prediction of basic properties of the SFG-population that 
relies on basic observables (e.g. the evolution of sSFR 
in main-sequence galaxies or their stellar mass distribu- 
tion) and uses their mathematical description (e.g. the 
Schechter function parametrization of the stellar mass 
distribution or slope and normalization of the main se- 
quence) to produce an analytico-empirical description of 
the statistical properties of SFGs. It can be both predic- 
tive (see, e.g., the indirect measurement of the evolution 
of molecular gas mass functions in Paper II) or help to 
(re-) interpret existing measurement s (e.g. IR luminos- 
ity functions or source counts; see iSareent et al.l 120121 
iBethermin et al.ll2012l iGruppioni et al.ll2013| ). ~ 
We introduced the 2-SFM framework in S12 where 
we demonstrated that the observational constraints on 
the z<2.5 IR luminosity functions can be reproduced 
based on only three observables: (i) the redshift evo- 
lution of the stellar mass function for SFGs, (ii) the 
evolution of the specific star formation rate (sSFR) of 
main-sequence galaxies, and (iii) a double log-normal 
decomposition of the sSFR-distribution at fixed stel- 
lar mass into contributions (assumed redshift- and 
mass-invariant) from main-sequence and starburst ac- 
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tivity. The split into (overlapping but offset) (s)SFR- 
distributions associated with main-sequence and star- 
burst activity is based on the distributions of sSFR pub- 
lished for massive (M+fMp > 10 10 ) SFGs at z ~ 2 pub- 
lished bv iRodighiero et alJ (|2011f ). The assumption that 
this double log-normal decomposition of the (s)SFR- 
distribution is invariant with stellar mass and redshift 
leads to a good agreement with IR-observables (a mild 
decrease of the importance of the starburst-component 
by <50% between z < 1 and leads to additional small 
impr ovements ; see discussion in S12 and lBethermin et alJ 

mil). 

The distinction between "normal" SFGs (implicitly as- 
sumed to be growing their stellar mass in a secular 
mode on the star-forming main sequence) and starburst- 
ing galaxies is central to the 2-SFM framework and of 
particular interest since it yields observationally verifi- 
able predictions of the notion that star-formation is a 
bi-modal process at low and high redshift. In this vein 
iBethermin et al.l (|2012f) assigned a characteri stic (albeit 
redshift-dependent; see IMaedis et all I2012IJ) IR spec- 
tral energy distribution to main-sequence and starburst 
galaxies and showed that this simple approach is ca- 
pable of reproducing the IR/radio source counts (inch 
new Herschel counts) at 24 to 1100 fim and 1.4 GHz. 
Given the sensitivity of the source counts, this obser- 
vation evidences that the 2-SFM framework provides a 
valid description of the dust emission from SFGs out to 
at least z~4, (i.e. over 84% of the age of the Universe). 
The 2-SFM description of the IR-properties of SFGs has 
also provided testable predictions which were verified in 
recent work, e.g. the redshift distrib ution of SCUBA- 
2 450 fj,m sources (jGeach et alJ |2012[ ) and the redshift 
distribution of lensed 1.4 m m sources detecte d with the 
South Pole Telescope (see iWeiss et al.l 120131 : Fig. 9). 
The good agreement of the predictions with the latter 
measurement suggests that the basic ingredients of the 
2-SFM framework (e.g. the minor role of starbursts) re- 
main applicable out to z ~ 6. 

4.2. Boosting of main-sequence galaxies: mathematical 
description 

Encouraged by the successful reproduction of the IR 
properties of the SFG-population we now further develop 
the 2-SFM framework with the primary aim of using it 
for a predictive analysis of the molecular gas properties of 
SFGs at high redshift. In preparation for this we revisit 
the key ingredient of the 2-SFM framework - the dou- 
ble log- normal decomposition of the (s)SFR-distribution 
at fixed stellar mass. Analogously to SI 2 we write this 
(s)SFR-distribution as the sum of two log- normal distri- 
bution functions Q describing main-sequence (MS) and 
starburst (SB) galaxies, respectively: 



p(sSFR)| M . = Sms(sSFR) +£ sb (sSFR) 



(5) 



Note that with this parametrization "starburstiness" is 
not an all-or-nothing property, but that this Ansatz nat- 
urally leads to a continuous spectrum of burst-bearing 
sources ranging from those with strongly boosted star 
formation activity to others with only a mild enhance- 
ment. 



Both the burst-bearing and the normal galaxy pop- 
ulation are described by an amplitude Ax with units 
of [Mpc~ 3 dcx(sSFR)" 1 ], a dispersion erx (units: 
[dex(sSFR)]) and a mode (sSFR) x (Xe{MS, SB}). In 
particular, the MS-distribution has the functional form 



Sms(sSFR) 
A M s exp 

A M s exp 



[log(sSFR) - log((sSFR) M s) 

zo MS 



2a 2 

ZCT MS 



(6) 



if, for the sake of brevity, we introduce a specific star 
formation rate x = log(sSFR/(sSFR)Ms) that is normal- 
ized to the stellar mass- and redshift-dependent average 
sSFR of main-sequence galaxies, (sSFR)ms- Similarly, 
for starbursting sources, wc write 

g SB (sSFR) = 



A SB exp 



A SB exp 



[log(sSFR)-log((sSFR) SB )]' 
[log(sSFR) - {log((sSFR) MS ) + B SB }f 



2^b 



[x - B SB ] 2 
A SB exp [ —2 



2a 



SL! 



(7) 



Here £?sb is the offset between the peak position of the 
main sequence and starburst component of the sSFR- 
distribution. In our interpretation it represents the av- 
erage sSFR-enhancement - or boost - brought about by 
the burst-inducing process. In the following we will as- 
sume that there is a process (to be discussed in Sect. I6.1[) 
which boosts the star formation activity of an (s)SFR- 
dependent fraction of main-sequence galaxies with ini- 
tial sSFR-distribution Gms- &MS i s identical with the 
log-normal distribution in eq. [51 except for a higher nor- 
malization A M S — > ^MS = ^MS x I (Sus+Qss:x)fl 

(Here /(/ : x) stands for the integral of the function / 
over the range x £ ] — oo, oo[.) 

The "boost-function" describes the spectrum of pertur- 
bations that main-sequence galaxies suffer. It is effec- 
tively a convolution kernel which transfers galaxies from 
the main sequence to the starburst-distribution: 



GsB{x) = (g% ss *BlC)(x) 



G^ s (y) BK.{x - y) dy 



(8) 



The boost-function kernel (BJC) is obtained using the 
convolution theorem which links the Fourier transforms^ 



Gsb(1c) = g^ s (k) x BlC(k) 



(9) 



15 For the one-dimensional Fourier transform and its inverse we 
use the following convention 



/oo 
- oc 

1 f°° - 
/(*) = o / /( fc K 



dx 

kx dk 



4.2.1. The boost-function: basic properties 



With this definition the Fourier transforms of the log-normal dis- 
tributions for starbursting sources and the unperturbed main se- 
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and by then applying the inverse Fourier transform: 



BfC(x) 



1 

2^ 



g SB (fc) 



ikx 



dk 



(10) 



Eqs. |9]and [TU] will not always have an analytical solution. 
In the present case however, since both parent and resul- 
tant distribution - and Gsb, resp. - are log- normal, 
the boost-function kernel also has this functional form: 

BK(x) = Qbk(x) = Cbk cxp (- [x ~ { 4 BK]2 ) ■ (11) 



By explicitly solving eq. [51 the three free parameters of 
the boost-function in eq. [TT] are found to be: 



a 



BK 



^sbosb 



1 



with 



MS ' 



and 



(12) 



SB 



The shape of the boost-function with free parameters 
given by expressions [12] is shown as a solid red line in Fig. 
Q] and compared to distributions of SFR-enhancements 
reported for simulated and observed interacting galax- 
ies. 

It is natural to expect that the process which statis- 
tically/physically links starbursting and normal galax- 
ies is galaxy interactions and merging. In Sect. 16.11 
we discuss in detail whether theory or observations can 
provide supporting evidence for such a straightforward 
connection. To summarize, some properties of the 2- 
SFM boost-function are suggestively reminiscent of SFR- 
enhancements in observations and simulations. Other 
aspects do not conform to the expectations of what a re- 
alistic boost distribution should look like if it accounts 
for, e.g. minor interactions, interactions including pas- 
sive galaxies and the fact that an observational snapshot 
of the starburst population will catch different objects in 
different phases of the burst. This difference, however, 
could be entirely due to the impossibility of statistically 
distinguishing between normal and only weakly star- 
bursting galaxies in a direct fit to the sSFR-distribution; 
the conventional view that starbursts are often tied to 
galaxy interactions hence remains a viable scenario and 
we now consider a modification to the boost-function 
that we expect to apply for the idealized case that all 
starburst events are triggered by merging. 

4.2.2. Statistical correction for paired, ante-merger galaxies 

In our presentation of the 2-SFM boost-function we 
have so far skipped issues that would complicate an 
immediate interpretation that is based purely on the 
mathematical description of the problem. A first and 
strong simplifying assumption is that the boost-function 



quence become 



Gss(k) = e 



MS k ~\ 2 



and 
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Di Matteo et al. (200B): 
1 ^ local mergers 

gas-rich mergers 

Scudder et al. (3012): 

SDSS galaxy pairs 

boost'— function 



20 



50 



excess (s)SFR 



Fig. 4. — : Comparison of the 2-SFM boost-function (i.e. 
the spectrum of (s)SFR-perturbations that move nor- 
mal galaxies off the main sequence and into a starburst- 
state) with measured and simulated distributions of 
SFR-enhancements induced by galaxy interactions. Red 
(solid) - 'direct' boost-function; red (dashes) - boost- 
function including explicit correction for merger statis- 
tics (see Sect. I4.2.2|) . In the cross-hatched area the 
shape of the 2-SFM boost-function is not known accu- 
rately (see discussion in Sect. 16. ip . The distribution 
of SFR-enhancements in massive (M* > 10 10 Mq) mem- 
bers of z € [0.02, 0.15 ] SDSS galaxy pairs measured by 
I Scudder et all (|2012t blue line) does not include mergers 
that have already undergone final coalescence and hence 
represents a lower observational limit to the total local 
SFR-excess distribution caused by interactions. Black 
and light grey lines - average of peak boosts (relative 
to isolated galaxies) me asured in SPH and g rid-based 
N-body simulations (cf. iDi Matteo et al.l 120081) of local 
mergers and gas-rich mergers with gas fraction similar 
to those of z ~ 2 main-sequence galaxies. 

is mass- an d redshift-independent. However, as discussed 
in S12 and lBethermin et al.1 ()2012[ ). current observations 
so far are consistent with this hypothesis. The sSFR- 
distributions of main-sequence galaxies and starburst- 
ing sources are snapshots that provide no direct infor- 
mation on the time scale over which star-formation in 
main-sequence galaxies is boosted (and over which they, 
possibly having undergone a merger, later fall back onto 
the relation or drop below it). The relative redshift- 
independence of the boost-function that is suggested by 
observations implies that the flux of galaxies into and 
out of the starburst-componcnt of the double log-normal 
distribution also should not evolve strongly with red- 
shift. If the SFR-enhancements were always the result 
of galaxy merging, then this would motivate a modifi- 
cation to the boost-function. We regard each pair of 
merging galaxies as a single system that is composed of 
two main-sequence galaxies. Because the parent sSFR- 
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distribution CJms(sSFR) is by assumption independent of 
mass and symmetric, both galaxies involved are drawn 
from the same distribution and their their average SFR 
will follow a log-normal function that is centered on the 
same sSFR as the original distribution C/ms( s SFR) but 
narrowed by a factor \/2. For the boost-function this im- 
plies a broader distribution of SFR-cnhancements while 
the peak location of the boost-spectrum remains identi- 
cal (see also distribution plotted with a dashed red line 
in Fig. |5): 



°BK = \/ ^ SB 

(x)bk =B S b 

Asbosb 

L'BK : 



and 

such that 



(13) 



1 



Al 



MS a MS V27TCTBK 



In the following we will refer to this version of the boost- 
function as "boost-function including an explicit correc- 
tion for mergers" (as opposed to the boost-function de- 
scribed by eqs. and henceforth called: "direct boost- 
function"). The principle findings of this paper are valid 
irrespective of the choice of boost-function, but for the 
sake of legibility we will only show results obtained with 
the direct boost-function where plotting both alterna- 
tives would reduce rather than improve clarity. 

4.2.3. The continuously varying boost-distribution 

As previously mentioned, the "typical" sSFR increase 
Bsb of starbursting sources is approx. a factor of four, 
but the average boost varies as a function of (s)SFR, as 
does the (relative) number of sources undergoing burst- 
like activity. For example, a source with measured 
sSFR twice as large as the (redshift- and stellar mass- 
dependent) main-sequence average could either display 
this excess simply due to a larger than average gas frac- 
tion and without having suffered any triggering, it could 
have experienced a modest boost, or - with a lesser prob- 
ability - it could initially have been a gas-poor, low- 
(s)SFR outlier to the main sequence which has been 
strongly boosted. In the following we quantify these vari- 
ations that are a consequence of the convolution in eq. 

m 

To find the "typical" boost of burst-bearing sources at a 
given sSFR (or log(sSFR /(sSFR)Ms) = xo) we consider 
the integrand in eq. 



Gms( x )Gbk(xq - x) 



'Ms( x 



&sSfr) £bk(&ssfr) 



(14) 

Here we introduced a variable for the logarithmic boost, 
frsSFR = £o — x = log(sSFR / sSFR) , in order to be able to 
directly locate the peak of the boost distribution, &"gp^, 
by solving the minimization problem: 



d 



db 



'sSFR 



Gms( x o - frsSFR.) £bk(& s sfr) \ = . (15) 



Given the properties of the exponential function this is 
equivalent to requiring 



d 



db 



sSFR 



(xq — bsSFB.Y 



2(T 2 
zo MS 



db s 



SFR 



(&sSFR — BqbY 

2^ 



= , 



(16) 



an equation which has the solution 



SB 



J sSFR 



1 + (" BK /»Ms)' 



(17) 



To determine the shape of the boost-spectrum which 
peaks at fo^FR we consider the product of the two Gaus- 
sians in eq. [M] 



Sms( x o — &ssfr) £bk(&sSfr) = 



(18) 



-4ms Cbk exp 



[xo 



sSFRj 



SFR 



BsbY 



2<s 



24k 



and examine the exponent 

[xp - fr s sFR] 2 [hsFR ~ -Bsb] 2 
2a 2 

which can be re-written as 



2^Ik 



r> (o'MSO'Bk) 2 
Z — ^ — T — 1 — 



^sSFR 



Xpg| K + ^SB^ms , 



sSFR 



'MS 



^O^BK 



'BK 



Bi 



SB^MS 



'MS 



(19) 



BK 



This is again a quadratic form, implying that the boosts 
at fixed normalized specific star formation rate xq are 

distributed log-normally with a width a^— ' (°" MSO " BK ) 2 



MS T BK 



that is independent of sSFR. 
On the left-hand side of Fig. [S] we visualize with differ- 
ent colors the changes in the boost-distribution for sS- 
FRs ranging from ~0.1 to two hundred times the char- 
acteristic main-sequence value, (sSFR)ms- Note that the 
sSFR-variations take place within a given bin of stellar 
mass and that the integral over the boost-distributions 
at all sSFRs would give a total boost-distribution that 
is equal to the boost-function plotted with the red solid 
line in Fig. |5 This is not immediately obvious based on 
Fig. [5] where we have scaled all sSFR-dependent boost- 
distributions such that they give the fraction of sources 
with boost 6 S SFR in a specific sSFR-bin. This repre- 
sentation highlights the evolution of the fraction f 8 of 
starbursting sources (given by the ratio of the two log- 
normal curves in eq. [5] see also Fig. (SJ right) with sSFR 
while simultaneously compensating for the variation of 
the total number of sources across the width of the main 
sequence. The amplitude of the log-normal boost distri- 
butions thus grows until sSFR/(sSFR)ms — 8 where the 
number of main-sequence galaxies becomes insignificant 
with respect to the number of starbursting sources (see 
also flattening of the evolution of f SB in the right-hand 
panel of Fig. [5]). 

An alternative quantity which traces the increas- 
ing importance of starburst-activity at successively 
higher sSFRs is the typical fractional contribution, 
CSB = (SFR-SFR M s, i0 it.)/sFR=l - KT haSFR , of the burst- 
induced activity to the total SFR of boosted sources. 
(Here SFRMs.init. is the SFR of the galaxy in the main- 
sequence state prior to boosting.) C is complementary 
to the starbursting fraction of the population, f 56 , in 
that it describes the impact of the boost on an individ- 
ual galaxy, while P B provides the number of galaxies of 
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Fig. 5. — : Illustration of the properties of the continuously varying, sSFR-dependent 2-SFM boost-distribution (i.e. 
the distribution of (s)SFR-cnhancements of starbursting galaxies; see Sect. 14.2.2)1 . Left - log-normal distribution of 
boosts (see eqs. H7lto[T9)) of SFGs with a specific sSFR. All distributions are normalized to the total number of sources 
at a given sSFR and colors vary according to the value of sSFR, where sSFR is referenced to the main-sequence average. 
For example, at sSFR/(sSFR)ms = 10 the median boost is approx. seven-fold and 95% of the starburst population have 
boosts in the range 2.5-17. Starbursts with sSFR/ (sSFR)ms = V 2 i on the other hand, have on average experienced weak 
boosting by only approx. 30% and 95% of such systems have boosts that are smaller than a factor three. The relative 
number of secularly evolving galaxies (plotted at boost = 1 with a bar of appropriate height) and of boosted sources 
is given by the burst-bearing fraction f 35 . Right - variation of the starbursting fraction f 35 with normalized sSFR, 
compared to the evolution of the typical fractional contribution, C SB , of the burst-induced activity to the total SFR of 
boosted sources at a given sSFR. (Here 'typical fractional contribution' is defined as the contribution of a source located 
at the sSFR-dependent peak of the boost-distribution.) Directly on the main-sequence locus (sSFR/(sSFR)ms = 1) 
starbursting sources are rare (f SB < 2%o) but in those systems which have experienced boosting the contribution of 
the burst-activity to the total SFR is significant (C SB ~ 47%). At an sSFR-excess of sSFR/ (sSFR)ms — 3, above which 
starbursts are more numerous than normal galaxies (f 56 > 0.5), the starburst-contribution to the total SFR of boosted 
sources is already clearly dominant (C SB ~ 71%). These numbers are based on the 'direct' boost-function that uses the 
best-fit parameters of the double log-normal sSFR-dccomposition in S12 (see also eq. [5] in this paper). 



a given sSFR within the 2-SFM framework that are sub- 
ject to such boosting. The relative variation of these two 
quantities is compared in the right-hand panel of FigJSJ 
Burst-bearing sources at the lower envelope of the 
main sequence are rare and typically have small SFR- 
enhancements (e.g. by 15% at -2 cms)- At an identical 
positive offset from the main-sequence locus (+2 cms) 
the star-formation activity of starbursting sources is on 
average triplecC3 with respect to their main-sequence 
state prior to boosting. Note that the lower tail of 
the 2-SFM boost-function formally allows for "negative" 
boosts (i.e. suppression of star formation activity in in- 
teractions) but this only occurs for a negligibly small 
fraction of <0.1%o and 2%o of the SFG-population for 
the direct and merger-corrected boost-function, respec- 

16 These numbers are for the boost-function including the ex- 
plicit correction for merger statistics. For the direct boost-function 
the typical boost fr^gp^ varies more slowly across the main sequence 
(see eq. [13 and Fig. [5J. 



tively. 

5. RESULTS 

In the following we will show how a simpler under- 
standing of the molecular gas properties of SFGs at all 
rcdshifts z < 3 emerges when variations of, e.g. SFE or 
gas fraction, about the typical value of main-sequence 
galaxies are considered. To be able to establish such 
normalized trends requires a reliable prescription for 
the evolution of slope and normalization of the star- 
forming main sequence with redshift. In Appendix [3] 
we parametrize the evolution of sSFR for main-sequence 
galaxies as a smoothly varying function of redshift and 
stellar mass (see eq. IA1|) which we fit to a compilation of 
sSFR-data from the literature. We find that on average 
a sSFR vs. M+ relation sSFRocMJ' with exponent 
0.2 reproduces the systematic shift between the sSFR- 
e volution of galaxies in different mass bins out to z ~ 3. 
The extrapolation of the sSFR-evolution to higher red- 
shift (as briefly proposed for the discussion of gas fraction 
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Fig. 6. — : Dependence of SFE (see scale on right) and gas depletion time Tdep. (scale on left) on SFR for individual 
galaxies in our calibration sample (panel a) and for stacked galaxies from IMagdis e t al.l (|2012bt panel b) . (All data 
and symbols as in Fig. [TJb and Ob). 



evolution in Sect. 15.2. 2p is speculative, since constraints 
on the shape of the main sequence at z > 3 are much 
sparser. 

We begin this section with our new description of SFE- 
variations between normal and starbursting galaxies 
(Sect. 15. ip and then discuss gas fractions variations and 
evolution plus empirical recipes for the CO-to-H2 con- 
version factor (Sects. 15.21 and 15.31 resp.). 

5.1. Simple recApes for star formation efficiency in 
massive star-forming galaxies 

5.1.1. Star formation efficiency in normal and 
star-bursting galaxies: observations 

The non-linearity of the integrated S-K law found in 
Sect. 13.2.21 implies a residual dependence of the gas 
depletion time, Td cp . = M mo i./SFR, and its inverse, the 
SFE, on SFR. Using our fit for main-sequence galaxies 
from cq. [3] we obtain 



= (a 2 , sfr-9) + log 



SFR 



o/yr 



M 



(20) 



= 0.18(±0.02) - 0.17(±0.03)xlog 



SFR 

M o/yr 



and 



log 



SFE 
Gyr-' 1 



(l-/3 2 )lo! 



/ SFR 

V M e/yr 



(0(2, SFR-9) (21) 



= 0.17(±0.03)xlog 



SFR 

M / yr 



0.18(±0.02) 



The dispersion about this characteristic value is ^0.2 dcx 
(see Sect. 13. 2j) . Fig. [5] illustrates how the galax- 
ies from our refer ence sample ( s ee Sect . I2.1[) and the 
stacked samples of IMagdis et al.l (|2012bl ) scatter around 
this average trend which - due to the general redshift- 



evolution of sSFR in SFGs - implies a roughly two- 
fold decrease of Tdep. between z ~ and z ~ 2 for galax- 
ies of M^/Mq =4xl0 4 which contribute most to the 



cosmic SFRD over this period (ICowie fc Bargerl 120081 
IGilbank et ail 120111 IKarim et al.l 120 111) . The variation 
between the depletion times of 1-2 Gyr in l ocal spiral 
galaxies (e.g. iLerov et al.l 120081 iBigiel et al.l 120111) and 
the approx. 0.5-1 Gyr determined for BM/BX- and BzK- 
selected galaxies at 1.5 < z < 2.5 (e.g. [Daddi et al.ir2010bl 
iTacconi et al.l [20Toh is much smaller than the difference 
between normal galaxies and strong starbursts (see off- 
set, dashed locus in Fig. [5]). 

We attempt to correct for the implicit redshift- 
dependence of SFE by considering a re-normalized effi- 
ciency. For each galaxy in our reference sample the nor- 
malization constant, (SFE)ms, is the SFE that a galaxy 
of equal gas mass would have if it lay directly on the 
inverse S-K relation given by eq. [3] In Fig. [7] we then 
plot the normalized SFE/ (SFE) ms a s a function of the 
normalized specific star formation rate sSFR/ (sSFR)mSi 
which in the 2-SFM framework is a good measure of 
starburstincss. The stellar mass- and redshift-dependent 
main-sequence average (sSFR)ms is calculated according 
to cq. [All 

With this choice for the representation of the data, main- 
sequence galaxies occupy the same region of parame- 
ter space regardless of their rcdshift. Our small refer- 
ence sample of starburst galaxies, on the other hand, 
is clearly offset from the main-sequence population in 
the plane of normalized SFE and sSFR. Note that we 
assign statistically-estimated CO-to-H2 conversion fac- 
tors based on a M*- and SFR-dependent metallicity 
Z(SFR, M*) to the majority of the normal galaxies in 
our reference sample, but that aco h as been directly 
measured for our sub-sample of starbursting galaxies. 
While the scatter of the normal galaxies about the main- 
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Fig. 7. — : Star- formation effici ency (SFE) vs . sSFR for selected main-sequence galaxies and starbursts at z<3 (left) 
and for stacked galaxies from IMagdis et al.l (|2012bt right). When normalized to the characteristic main-sequence 
value, (.)ms, a homogeneous behaviour of main-sequence galaxies at all redshifts becomes visible: SFEs vary very 
little within the main sequence while starbursting sources display enhanced SFEs that lead to their excess (s)SFR. 
(All data and symbols as in Fig. [5] In panel (b) we also show sample SFR-averages for local COLD GASS galaxies 
binned by sSFR-excess.) 



sequence average (SFE)ms is thus mo del- dependent, the 
strong SFE-excess found for starbursts is thus not an 
artifact of, e.g., assuming Milky Way- and ULIRG-likc 
conversion factors for main-sequence galaxies and star- 
bursts, respectively. Fig. [7J shows that the sSFR- and 
SFE-excess of the strongest starbursts in the sample are 
of a similar order of magnitude. This suggests that there 
is some kind of link between the SFR-cnhancement (or 
'boost' in the terminology of Sect. [4} and the increased 
SFE that starbursts display. The empirical calibration 
of this relation is the topic of the next section. 

5.1.2. Star formation efficiency in normal and 
star-bursting galaxies: the 2-SFM description 

Stellar mass and star formation rate are fundamental 
parameters in the 2-SFM description of SFGs. The tight 
and apparently redshift-independent integrated Schmidt- 
Kennicutt law found in Sect. 13.21 links the SFR and 
gas mass of main-sequence galaxies and hence provides 
a straightforward recipe to extend the 2-SFM framework 
to their molecular gas properties. At fixed stellar mass, 
in which case SFR/(SFR) M s = sSFR/(sSFR) M s we can 
write 

, / SFE \ , . f sSFR \ , . 

for the relation between normalized SFEs and sSFRs. 
This slow variation across the spread of the main se- 
quence with SFE cx sSFR°' 17±0 03 is superimposed on the 
data in Fig. [7b. 

The 2-SFM framework distinguishes between main- 
sequence galaxies and starbursting systems which sup- 
port an elevated level of star-formation activity com- 
pared to what is assumed to be an initial, pre-burst state 



where such galaxies were indistinguishable from the large 
population of secularly-evolving, normal SFGs. Hav- 
ing derived a prescription which links the SFE of main- 
sequence galaxies to their offset from the main sequence 
in eq. [55]we now seek a similar relation for the starburst- 
ing fraction of the population. We adopt the following 
paramctrization to describe the SFE of starbursts: 

log (sfSz) = ^ x 6sSFR (23) 

where 6 s sfr is the logarithmic boost introduced in eq. 1141 
and SFEMS,init. is the SFE in the main-sequence state, 
prior to the onset of the burst-activity. Since we re- 
fer the SFE to this initial state by definition no ad- 
ditional normalization constant is required in eq. 1231 
Observationally, the amount of boosting that the star- 
burst galaxies in our reference sample have experienced 
is obviously unknown. We shall thus assume that their 
SFR-enhancements correspond to the median boost (i.e. 
the peak location of the boost-distribution, 6^gf? R ; see 
eq. [T7|) which is expected for sources with an sSFR- 
excess as determined for these starbursts. In Sect. 14.2.31 
we derived the boost-spectrum at fixed sSFR-excess (or 
deficit), sSFR/(sSFR)mSj and calculated the shifting of 
its peak &"gp x R with the normalized sSFR. We reproduce 
the average trends in Fig. [5£i for both the direct boost- 
function and the boost-function including an explicit cor- 
rection for mergers. Since the latter scenario assumes an 
sSFR-distribution of paired, ante-merger galaxies that is 
narrower (see explanations in Sect. 14.2. 2[) . a larger boost 
is required on average to reach a given sSFR-excess. This 
fact is reflected in a steeper slope of the corresponding 
boost vs. sSFR/(sSFR)ms relation in Fig. [8}i. The 
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Fig. 8. — : Illustration of the two steps needed to establish a link between the star formation efficiency SFE of 
starbursting galaxies (normalized to the average value of main-sequence galaxies with a given H2-mass, (SFE)ms) and 
their (s)SFR-boost within the 2-SFM framework, (a) Variation of the 'typical boost' (defined as the peak of the sSFR- 
dependent spectra of boosts shown in Fig. [5]) of the burst-bearing fraction of SFGs with their sSFR (normalized to the 
M+- and redshift-dependent value of main-sequence galaxies, (sSFR)ms)- Two cases are considered: (1) boost-function 
including explicit correction for merger statistics {grey shading, cf. Sect. I4.2.2j) . and (2) 'direct' boost-function {hatched 
area). The shaded/hatched areas mark the 1 cr-error on the typical boost, and were derived based on the uncertainties 
associated with the decomposition of the sSFR-distributions of massive z ~ 2 SFGs into two log-normal components 
(for main-sequence and starbursting galaxies, resp.; see S12 or the schematic representation in Fig. [S]). Bold lines at 
the core of the highlighted confidence regions trace the best-fit variation of the typical boost values. 
(6) Empirical calibration of the relation between SFE-enhancement and boost amplitude u sing a reference sample of 
selected starburst galaxies with measured aco (|Solomon et al.l 119971 iMagdis et al.l l2012bf ) and the relation between 
normalized sSFR and the typical boost of panel (a). The reference sample of starbursts is located with respect to 
the x-axis using the trend lines of panel (a) for the 2-SFM boost- function including explicit correction for merger 
statistics {dark grey symbols ) and for the 'direct' 2-SF M boost- function {bla ck symbols). (Crosses/dots are used for 
local/high- z starbursts from iSolomon et alj (|1997|) and Magdis et al.l (|2012bf ). respectively; lsigma-error bars plotted 
account for observational uncertainty on SFR, M+ and M mo \. but not systematic uncertainties, see discussion in Sect. 
15.1.21 ) The shaded/hatched areas span the 68% confidence region for a power-law relation between excess SFE and 
boost amplitude as parametrized in eq. [23] and passing through the main-sequence locus highlighted schematically. 



shaded/hatched areas in this figure panel indicate the 
uncertainty on the average relation between boost and 
normalized sSFR, estimated with a full accounting of the 
errors on (and covariance between) the parameters of the 
2-SFM double log-normal decomposition of the sSFR- 
distribution (see also Fig. 1 in S12). 
In Fig. |5Jd wc show the result of assigning representative 
boost values to the starbursts in our sample and then 
fitting eq. 1231 In accordance with our assumption that 
they experienced the median burst expected for an ob- 
ject with their sSFR-cxcess, we equate SFE/(SFE)ms ~ 
the SFE normalized to the average main-sequence value 
- and SFE/SFE M s,init. - the SFE-excess with respect 
to the pre-boost, main-sequence state of each individual 
galaxy - for these sources. The j/-axis values of the star- 
bursts in Fig. |8p are thus identical to those in Fig. [7] We 
emphasize that our small reference sample of starbursts 
does not provide sufficient statistics to justify the func- 
tional form of eq. [23] Here we simply use this data to 
derive the best fit given the preceding choice of a plau- 



sible paramctrization; in this context, the linear rela- 
tion proposed in eq. [23] is the simplest possible form 
that can be envisaged. The best-fitting values of the 
slope are 7sfe = l-63+ ; 14 and 1.50+ ; 10 for the direct 
and merger-corrected boost-function, respectively. The 
quoted 1 a errors reflect the observational uncertainty 
on SFR, M* and M mo \. (i.e. L' co & «co) ~ as plot- 
ted in Fig. [5}d for our reference starbursts - but not 
the systematic uncertainties related to the calibration of 
the average sSFR and SFE of main-sequence galaxies, 
nor those related to the f unctional form of eq. [2 3l It 
is interesting to note that iDi Matteo et al.l (|2007t ) find 
SFE-enhancements in merger simulations that exceed the 
boost in SFR. The behavior of the simulations thus qual- 
itatively matches the supra-linear relation inferred here, 
according to which SFE/SFE M s,mit. oc (boost) 7SFE with 
7sfe > 1- 

A non-linearity of this kind is virtually inevitable in order 
to self-consistently match the SFEs of strong starbursts 
at the sSFR they display. If, instead, we were to describe 
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Fig. 9. — : Schematic comparison of the sSFR-distribution at fixed stellar mass (a) with the SFE-distribution at fixed 
molecular gas mass (b). The curves are representative examples, computed using the best- fit parameters of the double 
log-normal decomposition of the sSFR-distribution (eq. [5]) and of the power-law relation fitted to the dependence of 
SFE on boost-amplitude for starbursting galaxies (cq. I2li]) . Associated uncertainties arc not shown. 



the SFE of starbursts as a linear mixing of two distinct 
S-K laws (one for normal galaxies and one for starbursts) 
based on the starburst-contribution to the total SFR of 
boosted sources (C SB , see Fig. [5p), the high SFEs ob- 
served for, e.g., local ULIRGs could only be reproduced 
when assuming a 'template' S-K law for starbursts which 
is offset to higher SFE than observationally seen for even 
the most extreme sources. 

An immediate consequence of the inferred supra-linear 
relation between SFE- and SFR-enhanccments in star- 
burst events is that the distributions of the SFG popu- 
lation with respect to (s)SFR and SFE are qualitatively 
different. (Here we consider the SFE-distribution at fixed 
gas mass.) We illustrate this in Fig. [9] for the case of 
the best-fitting double log-normal decomposition (eqs. [5] 
and [7]) and SFE-boost relation for starbursts (eq. [23]) . 
While the sSFR-distributions of main-sequence galaxies 
and starbursts blend, the two sub-populations are more 
strongly separated in terms of SFE. In comparison with 
the distribution of galaxies in the M mo i. vs. SFR plane 
in Fig. [T|d, we see that the bimodality predicted by 
the 2-SFM framework is less pronounced. So far, CO 
follow-up observations of SFGs have g enerally explic- 
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)lll iMagdis et all l2012al) . It is thus very likely that 
this selective observing strategy has artificially deepened 
the expected trough between the main-sequence and star- 
burst component in Fig. [9Jd into the genuine gap that is 
seen in Fig. QJ). 

We now have all ingredients to provide an empirical pre- 
scription for the SFE of starbursting systems. In analogy 
to the expression for normal galaxies in eq. 1221 we write 



the recipe for starbursts in terms of a normalized SFE, 

/ SFE \ / SFE \ /SFE MS ,init. 
lo S TFE^v = log ^ + lo § /orm 1 



V (SFE) MS 

= 7sfe x 6 sSFR + (1-/Ja) x log 



SFEMS.init. / \ (SFE)mS 

sSFR ms . 

init. 



(sSFR) 



MS 



(24) 



which is a combination of the SFE of the initial, main- 
sequence state prior to boosting (second summand; see 
also eq. l22j) and the burst-induced enhancement of this 
initial SFE (first summand; see also eq. |23|) . It is im- 
portant to realize that eq. [Ml stands for SFE-changes in 
individual starbursts rather than describing the sSFR- 
dependence of the average SFE-excess of the whole star- 
bursting population. The latter trend will be discussed 
in the next paragraph. 

By incorporating the prescriptions in eq. [22] and [24] we 
can predict the variation of SFE throughout the entire 
M* vs. SFR plane and for the entire SFG population (as 
opposed to individually for the sub-population of nor- 
mal galaxies and starbursts). In doing so, we will again 
make the simplifying assumption that all recipes are in- 
dependent of stellar mass, which reduces the problem to 
a calculation of the evolution of SFE with (s)SFR. The 
prediction involves several steps which we detail here in 
bullctizcd format for maximal clarity and in preparation 
of analogous procedures for the variation of gas fractions 
and aco m Sects. 15.21 and 15.31 respectively: 



• The hypothetical population of SFGs is distributed 
in M* and SFR according to the double log-normal 
distribution in eq. [5] 

• At each point in the M* vs. SFR plane an sSFR- 
dependent fraction f 50 of galaxies will fall into the 



18 



Sargent et al. 



100 



H 10 

fa 
m 



0.1 7 



■ main sequence 
. . • • ■ starbursts 
— all SFGs 




_ SFE « sSFR"' 



(a) 



10 



i 
o 

J 



0.1 



- 1 — I — I I I I I 



« sSFR" 




(b) 



0.1 



10 



0.1 



10 



sSFR/(sSFR) u 



sSFR/(sSFR) u 



Fig. 10. — : Average variation of normalized SFE (left) and molecular gas mass fraction /i mo i. =M mo i./M* (right) with 
normalized sSFR, as predicted by the 2-SFM framework. Data points plotted in grey in the background are as in 
Figs. [7] and [TT] Dashed line - expected average trend for main-sequence galaxies (sSFR-dependence as annotated 
adjacent to line); dotted line - expected average trend for the starbursting fraction of the star- forming population 
('direct' boost-function assumed); solid line - average evolution for the total population of SFGs; light grey shading - 
expected 1 a scatter around average SFE- variation. The average evolution of the total population traces the variation 
of the median SFE (/i mo i.) of the combined SFE-distribution (/i mo i. -distribution) of normal and starbursting galaxies. 
Its variation thus reflects the relative importance of main-sequence and starburst galaxies with offset from the main- 
sequence locus (sSFR)ms- When normalized to the typical main-sequence value, the predicted trends do not depend 
on redshift due to the simple power-law relations between SFR and molecular gas mass that are assumed in the 2-SFM 
framework (see Sects. I5.1.2l and l5.2.1l for details). 



starburst-category (see Fig. [51 right) and hence 
require a different recipe for the computation of 
SFE than is applied to main-sequence sources. 

• For main-sequence galaxies of a given sSFR, the 
SFE is computed with eq. [22] 

• At each sSFR, relations [T7] to [H from Sect. ET231 
allow us to construct the spectrum of boosts 
&sSFR and to infer the former pre-burst efficiencies, 
SFEMS,init., of the burst-bearing systems. Their 
SFEs then follow from eq. El 

• A 'typical' SFE - here we use the median - is calcu- 
lated for the joint population of normal and star- 
bursting galaxies. It reflects the relative impor- 
tance of the two sub-populations at a given location 
in the M* vs. SFR plane. 

Fig. [TUk shows how, beginning at the lower edge of the 
main-sequence locus, the median SFE of the total SFG 
population first traces the slow rise of SFE which char- 
acterizes normal galaxies. Meanwhile, the SFE of burst- 
bearing sources (dotted line) rises more quickl^vFI but 
such starbursts are exceedingly rare within most of the 

17 The average SFE of starbursting systems is predicted to 
have a slightly sub-linear dependence on sSFR. This is the conse- 
quence of convolvin g the supra-linear evolution of SFE with boost- 
amplitude (see eq. I23J I with the shallow dependence of average 
boost amplitude on sSFR (see Fig. [SJ. 



main sequence. They become the dominant component 
of the SFG population at around 3 < sSFR/(sSFR) M s < 4 
where the global median abruptly jumps to join the trend 
of steadily rising SFE of starbursts, which reflects the in- 
creasing boost-amplitudes that are required to reach the 
highest sSFRs. Even though the 2-SFM framework as- 
sumes a full continuum of SFR- and SFE-enhancements, 
the changing population mix between normal and star- 
bursting galaxies thus has lead to a "bimodal" behaviour 
of the SFE. The exact shape of the jump depends not 
only on the parameters of the double log-normal decom- 
position but also on the SFE-dispersion of main-sequence 
galaxies. In plotting the average trend in Fig. [TCTk we 
have assumed a scatter of 0.2 dex, in accordance with the 
measurement made on Fig. Q] The scatter (indicated 
with light grey shading in Fig. I10[) abruptly increases in 
the transition region with its mixture of main sequence 
and starburst galaxies and then at sSFR/(sSFR)ms ^ 4 
is predicted to settle to a constant value that is somewhat 
larger than the 0.2 dex of the main-sequence locus as it 
simultaneously reflects (a) the spectrum of initial main- 
sequence states that end up at a given sSFR-excess by 
virtue of their different SFR-boosts, and (b) of the dis- 
persion in SFEs in the pre-burst main-sequence state. 
The depicted average trend also depends on the relation 
between SFE and boost. In Fig. [TOa we show the predic- 
tion for the case of the 'direct' boost-function (the SFE- 
evolution for starbursting galaxies would steepen when 
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Fig. 11. — : Molecular gas mass fraction (/x mo i. = M mo \_ /M+) vs. sSFR f or selected main-sequence galaxies and star- 
bursts at z<3 (left) and for stacked galaxies from lMagdis et all (|2012bl : right) (all measurements normalized to the 
M±- and redshift-dependent average of the main-sequence and all symbols and data as in Fig. [6]). The gas fractions 
of normal galaxies rise uniformly (independent of redshift) across the main sequence while starbursting sources have 
a gas content that is somewhat lower than the main-sequence average. 



using the boost-function including an explicit correction 
for mergers) and do so only for the best-fit value of 7sfe- 
We refrain from estimating formal errors for this aver- 
age trend since the systematic uncertainties (for exam- 
ple those pertaining to the choice of a function relating 
SFE-excess and boost as in eq. 123]) in any case strongly 
outweigh these. 

Note that, by working with (s)SFRs that are normalized 
to the average main-sequence value rather than absolute 
quantities, we have removed all dependence on redshift. 
Our prediction for the variation of SFE can thus be sum- 
marized by a single track, as shown in Fig. 110a . It is gen- 
eralizable to the full range of redshifts and stellar masses 
by (1) multiplying (s)sSFRs by the normalization con- 
stant, (sSFR)ms, which is given by cq. IA11 and by (2) 
using eq. [22] with SFR = ill* x (sSFR) M s to obtain the 
normalization, (SFE)ms- an d thence the absolute value 
of SFE. 

5.2. Molecular gas fractions in star-forming galaxies 

5.2.1. Gas fractions in normal and star-bursting galaxies 

The predictions of the previous section for the av- 
erage variation of SFE in SFGs with offset from the 
main-sequence locus are equivalent to a variation of the 

molecular gas mass M mo i. = {M* x sSFR(M*,z)}/SFE. 

By construction, the stellar mass M* in this expres- 
sion is known in the 2-SFM approach, implying that we 
can directly compute the molecular gas mass to stellar 
mass ratio, ^ mo i. = Mmd./M*, as a function of sSFR. 
To obtain predictions that are independent of redshift 
and stellar mass we again consider normalized quantities, 
sSFR/(sSFR) MS and /w. /(Mmoi.), in Fig. [TUb. The bi- 
modal behaviour of the average SFE-evolution leads to 
two distinct regimes: (1) a nearly linear increase of the 



gas fraction across the main sequence, and (2) for strong 
starbursts, gas fractions that vary more slowly with sSFR 
and are somewhat lower than the gas fraction of the av- 
erage main-sequence galaxy. (Note that as in Fig. [POk . 
gas fraction variations within the starbursting popula- 
tion are again shown for the direct boost-function; using 
the boost-function that includes the explicit correction 
for mergers would result in a shallower trend.) As in the 
case of SFE (see Fig. [TOk). the transition between these 
two regimes is almost step-like and characterized by a 
large dispersion in yU mo i.. We explore the link between 
gas fractions during a starburst episode and prior to the 
onset of burst-activity in more detail in section 16.31 
We compare these expectations of the 2-SFM framework 
with real data in Fig. [TT] While our reference sample 
of starbursts is too small to quantitatively constrain any 
residual variation of gas fractions at high sSFRs, the pre- 
dicted rise of gas fractions across the main sequence is 
well sampled by the reference sample of normal galaxies. 
As discussed in Sect. 15.2.21 the gas fractions inferred for 
90% of the main-sequence galaxies involve an assumption 
about the metallicity-dependence of aco- Nevertheless, 
they show no systematically different behaviour than the 
eight sources (boxed symbols in Fig. [TT]) for which mea- 
surements of ac o exist. A very simila r slope (sSFR 9 ) 
was measured bv lMagdis et al. (20129) for stacked sam- 
ples of main-sequence galaxies divided in bins of sSFR- 
excess in which the gas mass was constrained via the 
far-IR dust emission. Stacking-based measu rements of 
Mmoi. at z = l, and 2 bv lMagdis et al.l (|2012bf ) are shown 
in Fig. [Tib and found to coincide with the 2-SFM predic- 
tions and gas fractions determined on a n individual ba- 
sis for galaxies in our reference sample. ISaintonge et al.l 
were able to sample molecular gas mass fraction 
variations over a larger range in sSFR which extends 
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to significantly below the star forming main sequence 
and also slightly into the starburst-regimc. In their lo- 
cal COLD GASS data set the average /z mo i. scales as 
approx. sSFR 7 . A significant deviation from the 2- 
SFM predictions for normal galaxies is only seen in their 
highest sSFR-bin which lies in the transition region be- 
tween main sequence and starbursting outliers and may 
hence be expected to reflect the transition to the lower 
gas fractions in starbursting galaxies (see also the indi- 
cation of a SFE-increase in the COLD GASS data set at 
sSFR/(sSFR) MS =4.5 shown in Fig. Wp)- 

5.2.2. Gas fraction evolution across cosmic time 

The well-defined relations between M+ and SFR, and 
between SFR and H2-mass allow for a straightforward 
prediction of the redshift-evolution of the molecular gas 
fraction, / mo i., in normal galaxies. Using eq. |3]we can 
write 

Afmol. 1 



mol. : 



M mo i. + M* 1 + M*/(const.xSFRft») 
1 



f 



x sSFR-fe 



(25) 



where const. = 10 Q2 - SFR . If we insert the average 
sSFR(M 1t , z ) o f main-sequence galaxies (parametrized 
as in eq. IAI[) into this equation we obtain an evo- 
lutionary trend which we plot in the upper half of 
Fig. [12] for four different stellar masses in the range 
5x 10 9 < M*/M < 10 11 . The gas fractions of normal 
galaxies at z < 3 predicted in this manner are in ex- 
cellent agreement with literature data, which suggests 
that the gas content of secularly-evolving SFGs is an im- 
portant driver of the cosmic sSFR-evolution. Support- 
ing evidence for this tight link between the evolution 
of sSFR and the gas fraction of main-sequence galax- 
ies was recently provided by the analysis of more than 
50 SFGs at 1 < z < 3 w ith CO-fhix m e asure ments from 
the PHIBBS survey in iTacconi et all (120131 ; dark blue 
stars in our Fig. H2j). ICombes et all (12013ft . on the other 
hand, hold a joint redshift-evolution of both SFE and 
gas fractions by roughly similar amounts to be respon- 
sible for the cosmic evolution of star formation activity. 
This apparent discrepancy can be explained by the dis- 
tinctly different behaviour of normal galaxies and star- 
bursts with respect to SFE and / mo i., which we illustrate 
in Fig. [lOl a nd by the fact t hat the 0.2 < z< \ ULIRGs 
analyzed by ICombes et al.l (|201lL 12013ft have generall y 
large sSFR-excesses (see Fig. 9 in ICombes" etaHl2013h . 
while their local reference sample overlaps with the main- 
sequence population. These authors hence tend to com- 
pare fairly low-efficiency z = systems with z > 0.2 star- 
bursting galaxies (thereby overestimating the importance 
of SFE for sSFR-evolution) which on average have lower 
gas fractions than equally massive normal galaxies (lead- 
ing to an underestimate of the gas fraction evolution). 
The link between sSFR and gas fraction also manifests 
itself at fixed redshift as a variation of sSFR across the 
main sequence (see Fig. ITT]) . This is equivalent to stating 
that the dispersion of the main sequence can be at least 
partially ascribed to different gas fractions. In a recent 
morph ological study of z ~ 1 disk galaxies iSalmi et al.l 
(|2012ft reported that systems with clumpy substruc- 
ture are found to be systematically offset to higher val- 



ues of sSFR than their smoother counterparts. Since 
clumps arc a telltale signature of violent disk instabili- 
ties in gas-rich high-redshift galaxies (e.g. lAgertz et alj 
20091 iCeverino et al.ll2010llForster Schreiber et al.l l201lL 



Swinbank et al.l I2011L iWuvts et al.l 12012ft this observa- 



tion thus provides independent evidence for increasing 
gas fractions within the main sequence. 
Following the publication o f simulation-based recip es for 
the calculation of aco by iNaravanan et al.l (12011ft. the 
same authors have recently questioned (jNaravanan et aLl 
l2012f ) the reliability of the high gas fractions reported 
in the literature. (Alternative predictions for aco that 
make use of the observed relations between L' co and 



Lir and SFR and H 2 -mass are presented in Sect. I5.3H 
Various additional predictions for g as fraction evolu- 
tion based on nu merical simulations (jDave "eFall [2TTTTL 
Duffy et al.ll2012t ) or semi-analytical modeling ( Fu et all 
2012ft are shown in the lower half of Fig. [T2] Although 



these predictions generally lie within the range of gas 
fractions expected in the 2-SFM framework for galaxies 
located on the main sequence, a tendency for shallower 
redshift evolution of / mo i. in the simulations is ubiqui- 
tous. Rather than being the consequence of incompati- 
ble assumptions for the calculation of aco these differ- 
ences compared to the 2-SFM predictions might reflect 
the wel l-known problem that both semi-analytical mod- 
els (e.g. lFontanot et al.ll2009ft and cosmological h ydrody- 
namical simulations (e.g. iWeinmann et al.ll2012ft tend to 
produce too many stars too early in the history of the 
universe, especially in lower mass galaxies. The result- 
ing accelerated exhaustion of gas reservoirs would then 
likely lead to lower g as fractions than we predict using 
the 2-SFM approach. iPopping et"aTI (|2012ft used a simi- 
lar, empirically motivated approach as the one proposed 
here to indirectly infer the gas content of both late- and 
early- type galaxies at z < 2. Here we show the redshift 
evolution of / mo i. these authors derive for main-sequence 
galaxies with the stellar masses plotted in Fig. Q2J Their 
expectations are in good agreement with ours at z < 1 
but begin to differ from the 2-SFM predictions at higher 
redshift, probably due to incompleteness in their lower- 
mass (M^/Mq < 10 11 ) galaxy samples. 
We end this section by plotting explicitly in Fig. IT^lthc 
M*-dependence of the molecular gas fraction which was 
already visible in Fig. [12] as a vertical offset between the 
evolutionary trends. Two panels with the 2-SFM predic- 
tions for z ~ 1 and 2 are shown. In both cases the range 
of expected gas fractions / mo i. = M mo i./(M*+M mo i.) at 
fixed stellar mass can be significant, e.g. for a stellar 
mass of 5xlO lo M0 a A/ mo i. of approx. 0.5 is expected, 
depending on whether the source is located at a positive 
or negative offset of 2 oms with respect to the main se- 
quence. We illustrate this considerable scatter with our 
two statistical samples of GOODS-S galaxies (see Sect. 
12. 2p . to which we apply eq. [3] to indirectly infer gas 
masses. In spite of the large dispersion, a clear trend 
of decreasing gas fractions with increasing stellar mass 
is seen. Predictions from semi-analytical models show 
qualitatively similar, albeit somewhat shallower trends. 

5.3. The CO-to-H-i conversion factor aco 

In Sect. 15.1.21 we showed how the sSFR-dependcnt 
starburst demographics of the 2-SFM approach led to 
a nearly step-like variation of the SFE even if star- 
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Fig. 12. — : Redshift-evolution of the molecular gas fraction, / mo i., of main-sequence galaxies in four different stellar 
mass scales, as predicted by the 2-SFM framework based on the evolution of the sSFR and of the integrated S-K 
relation (see eq. I25[) . Upper half of figure - comparison with literature data (taken to have stellar masses within at 
most a factor two of the mass scale used for the analytical predictions) ; lower half - comparison with predictions from 
numerical simulations and semi-analytical modelling. T he d ark grey shading illustrates the uncertainty (±1 ct; reflects 
the uncertainty of the sSFR-evolution according to cq. IA1I and of the integrated S-K law in cq. [3]) on the evolution 
of /moi. for a typical main-sequence galaxy. Light grey areas illustrate the predicted dispersion of gas fractions. The 
evolution for galaxies offset by +2/+l/0/-l/-2<7 from the average main-sequence locus arc additionally highlighted 
(uppermost to lowermost medium grey line). The good agreement between predictions and data in the upper half of 
the figure - i.e. the observation of a synchronous evolution of sSFR and / mo i. - is consistent with the evolution of the 
gas reservoirs in normal galaxies being the primary driver of the cosmic sSFR evol ution. Star-shaped symbols indicate 
gas fr actions determined with stacking ([Magdis et al.ll2012bT ) or sample averaging ([Saintonge et al.|[2012l iTacconi et al.l 
120131) . Colors/line styles used to represent different simulation predictions in the lower half of the figure are: dotted/dot- 
dashed / dashed / so lid green lines for the vzw/cw/nw/sw (momentum-conserving / constant /no / slow w ind) scenario in 
iDave et al.l ()2011l ): red/scarlet lines for th e 'LOSONSI^Y'L IOONS^' realizations in lDuffv et ail (|2012h : light/dark blue 
lines for 'prescription l'/'prescription 2' in iFu et ahl (|2012t H2 fraction depending on local cold gas surface density and 
metallicity/H2 fraction depending on ISM pressure), and for the assumption of two different star-formation laws for 
regions where atomic and molecular gas dominates, resp. (see iBigiel et al.l [20081 ) . 
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Fig. 13. — : Variation of the molecular gas fraction, / mo i., 
of main-sequence galaxies with stellar mass in two dif- 
ferent redshift bins [z ~ 1 - top; z ~ 2 - bottom) . Se- 
lect ed main-sequence gala xies from iDaddi et al.l (|2010al ) 
and iTacconi et al.l (|2010t) . as well as ave rage gas frac- 
tions de termined in the stacking analysis of lMagdis et al.l 
(|2012bD are plotted in color. Predictions from the 2-SFM 
framework and recent theoretical/numerical work in the 
literature arc plotted with the identical symbols and color 
scheme as used in Fig. [121 Black dots - indirect mea- 
surements of /moi. for 24 /xm-selected GOODS-S gal axies 
based on the inverse S-K relation calibrated in Sect. 13.21 
The hatched line indicates the completeness limit of the 
statistical GOODS-S sample, as determined by the depth 
of the 24 (im imaging and the average redshift of the two 
bins displayed. 



bursting galaxies are treated as a continuous exten- 
sion of normal galaxies, with depletion times that de- 
crease in proportion to the burst strength (referred to 
as 'boost' throughout this paper). The actual value of 
SFE is tightly linked to the CO-to-H2 conversion fac- 
tor, qco- In the past it has been common practice to 
adopt one of two discrete, 'consensus' values when es- 
timating molecular gas masses for high-rcdshift galax- 
ies: Qco =4.4M Q (Kkm/spc 2 ) -1 (the conversion factor 
which is found to appl y to Giant Molecul a r Clouds in 
the M ilky Way; e.g. [Bolatto et all 120081 lAbdo et al.l 
120101 ) for normal galaxies and 0.8M Q (Kkm/spc 2 ) -1 , 
a rep resentative average for lo cal starbursting ULIRGs 
(e.g. iDownes fe Solomon] 119981 ). Although the collec- 
tively high luminosities of distant, CO-detected galaxies 
suggested that their aco-values should be ULIRG-likc, 
the first actual estimates of aco for z > 1 disk galaxies 
(based on both dynamical arguments as in IDaddi et al.l 



I2010al or on the gas-t o-dust ratio approach implemented 
in lMagdis et al J 120111) were all broa dly consistent with a 
Galactic aco- Going a step further. [Genzel et al.l ([2012D 
subsequently were able to show that the conversion factor 
of high-z main-sequence galaxies scales with gas-phase 
metallicity in a similar manner as the negative power 
laws observed for local galaxies (e.g. iWilson 119951 Hsrael 



2012) 



laws observed tor local galaxies (e.g. IWnsonllli^a llsrae 
i~997l iBoselli et al.l2002tlLerov et al.ll2011llSchruba et al 



The physics of the multiphase ISM which ultimately 
determines the exact value of aco is complicated, re- 
gardless of whether the emission from individ ual star- 
forming regions (e.g. iGlover fc Mac Low! 120111) or from 
larger sca les even up to integrated emission are consid- 
ered (e.g. iNaravanan et all 120111 iFeldmann et al.1120121 
iPapadopoulos et al.1 I2012T ) A common feature of all 
these theoretical or numerical calculations is a depen- 
dence of aco on metallicity. Here we adopt a scaling 
aco o c Z~ 9 for main-sequence galaxies, as motivated in 
Sect. 13.21 In the general spirit of the 2-SFM approach, 
the aco-recipe proposed here for starbursting sources 
does not attempt to capture the aforementioned complex 
ISM physics in its entirety, but rather treats it implic- 
itly by letting the data dictate the functional form of the 
recipe. The fact that starbursts in the 2-SFM framework 
preserve a memory of their former main-sequence state 
means that our predicted starburst-aco values also de- 
pend on metallicity but in a more complicated way which 
is detailed in the following. 

5.3.1. Conversion factors for starbursts: empirical 
calibration of boost- dependence 



As for the SFE (see eq. [22]) we assume that aco varies 
smoothljH with the boost of a starburst, 



log 



/ QCO 
\ a MS, 



init. 



7« C 



x b 



sSFR 



(26) 



and use our sample of reference-starbursts (cf. I2.1.3[) 
to determine the most suitable value of 7 QC o > given 
this choice of parametrization. Boosts 6 s sfr are as- 
signed as in Sect. 15.1.21 and aco-values are give n in 
IDownes fc Solomon! (|1998l ) or iMagdis et all (|2012bf ) for 
each starburst in the reference sample, oms, init. cor- 
responds to the conversion factor of a main-sequence 
galaxy with the same SFR as a reference-starburst, but 
with L'qq and A7h 2 given by the inverse S-K relations 
in eq. []] and [3] respectively. Since we refer the aco 
to this initial state by definition no constant term is re- 
quired in eq. Solving for 7q, co we obtain 7 QCO = - 



0.78 



+0.07 



(la errors quoted) for the case of the direct 



boost-function and 7 QCO = -0.72±0.05 for the merger- 
corrected boost-function. The boost-dependence for the 
two scenarios is shown in Fig. H"4"k . together with the 

18 The functional fo rm of eq. 1261 is not merely motivated by 
its symmetry with eq. 1231 it also reflects the expectation that 
the state of the ISM evolves continuously, e.g. depending on the 
strength of the tidal forces which may enhance the amplitude 
of the turbulent motio ns during galaxy-galaxy interactions (e.g. 
Bournaud et al. 2011b). Larger velocity gradients and higher tem- 
peratures, which are generally characteristic of the turbulent and 
dense starbursting ISM, both lower aco while higher column den- 
si ties increase the the CO-to -H2 conversion factor (see, e.g., eq. A4 
in Papadopoulos et al. 2012). 
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Fig. 14. — : (a) Empirical calibration of the decrement in aco with the amount of (s)SFR-boosting experienced by star- 
bursts (all symbols and data as in Fig. [5b). The decrement is referred to the CO-to-H 2 conversion factor, (o:co)ms, init.i 
that would be expected for each starbu rst galaxy if it were a secularly-evolving, average, normal galaxy with the same 
molecular ga s fraction (see Sect. [5.3.11 f or details). Starburst aco- m easurements plotted i n black are from dynamica l 
modelling in lDownes fc Solomon! (|199*8l ). light grey symbols are for aco-values derived bv IPapadopoulos et al.1 ( 20121 ) 
using a two-phase LV G model. (Crosses/dots are used for local/high- z starbursts from ISolomon et al I (|1997lT and 



iMagdis et al.l (|2012bl ). respectively; lsigma-error bars plotted.) The shaded/hatched areas span the 68% confidence 
region for a power- law relati on between starburst ago a nd boost amplitude as parametrized in cq. [2E] and when fitting 
to aco-measurements from lDownes fc Solomon! (119981 ). The light grey dashed line shows the best-fitting power- law 
relation inferred based on the aco-values from two-phase LVG modelling. 

(6) Predicted variation of average, normalized aco with sSFR-cxccss for normal main-sequence galaxies {dashed lines), 
starbursts {dotted lines] 'direct' boost-function assumed) and the total population of SFGs {solid lines). Data points 
plotted in grey in the background are as in Fig. 1101 Du e to the non-linear de pendence of the gas-phase metallicity 
on SFR and stellar mass (as parametrized by the FMR of lMannucci et aLll2010|) the average aco-trends predicted by 
the 2-SFM framework are both stellar mass- and redshift dependent. In a given redshift bin (color-coded as shown in 
the upper right corner) the shallowest variation across the main sequence (i.e. in the range 1 /& < sSFR/ (sSFR)ms < 6) 
occurs for the most massive of the three stellar mass bins considered here {M^/Mq = 5x 10 9 , 5xl0 10 and 5xlO n ), 
while the steepest variation occurs for the least massive bin. At high excesses of sSFR the aco-values are predicted 
to decrease the most steeply in the highest and most slowly in the lowest stellar mass bin plotted. Over the range 
< z < 2 shown here, the variation with stellar mass is expected to be more significant than with redshift. At the 
highest sSFR-excesses (boosts) the aco-values in starbursts plausibly asymptotically approach the lower limit set by 
optically thin CO line-emission. 



starburst-data used for the fit. The relatively slow, sub- 
linear decline of aco with boost-amplitude (i.e. SFR- 
cnhancement) according to aco (boost) 7a co implies 
that to reach ULIRG-like values of the conversion fac- 
tor, which are about 1 /b of the typically assumed Milky 
Way-like aco = 4.4 A/ (Kkm/spc 2 ) -1 , a boost by a fac- 
tor of eight to ten is expected according to the 2-SFM 
description. 

Several of the aco _m e asur cmcnts for our reference star- 
bursts deviate more strongly from the average trend be- 
tween boost and conversion factor than was the case 
for the relation between SFE and boost we calibrated 
in Fig. [Hb- As a consistency check we hence used 
the results of the LVG ra diative transfer modelling by 
IPapadopoulos et al.l (|2012D of all nine starbursting lo- 
cal ULIRGs in our reference sample to re-derive the 
logarithmic slope 7 QCO in eq. [26[ We find that 



aco-values determined with one-phase radiative trans- 
fer models ar e on average consist e nt wit h the dynamical 
estimates of iDownes fc Solomon! (|1998f ). such that the 
resulting slope is almost identical to the previously mea- 
sured one: 7 QCO =-0.81^Qgg (direct boost-function) and 
7«co =-0.75±0.06 (boost-function corrected for merger- 
statistics). CO-to-fi2 conversion factors inferred with 
two-phase (for high- and low-excitation gas) ISM models 
are generally higher (see light grey crosses in Fig. [T4k). 
leading to a shallower slope j aco = -0.48^;?^ (-0A4±%°£) 
for the direct (merger-corrected) boost-function. Given 
the good agreement between the former two estimates of 
7«co we have adopted the dynamic ally-constrained CO- 
toTT conversion factors reported in IDownes fc Solomon! 
(|1998| ) throughout this paper. 

We can now write an expression that relates the conver- 
sion factor of a starburst in general to the main-sequence 
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average: 

, / "CO \ , / «CO \ , , / "MS.init. \ 

log 7 r — = log + log >- — 

\\ a Co)MS/ VOMS.init./ \(aCO/MS/ 

= 7«co x b s SFR + — 

... + £x [/(SFR MS ,init.)-/((SFR) M s)] ■ (27) 

Here the first term - which describes the aco-deficit of 
starbursts with respect to the pre-boost, main-sequence 
state of each individual galaxy - corresponds to eq. 1261 
The second term relates the pre-boost, main-sequence 
state to the main-sequence average and depends on the 
slope £ of the relation between «co and mctallicity (eq. 
[2]). Eq. [27] strongly resembles eq. [24] for the normalized 
SFE of starbursts but has the subtle difference that the 
term describing variations within the main sequence, 

log ( , a °? j = f X log ( — ^— ) 

= £><[/ (MO.32) -/ ((M0.32)ms)] 

= Cx[/(SFR)-/((SFR) MS )] , (28) 

has higher order terms in log(SFR) due to its de- 
pendence on the FMR parameter /xn 32 = lp g(Af^)- 
0.32xlog(SFR) defined in IMannucci et all ([20101) . Since 
/(SFR) = C[log(SFR) 4 ] in general (see Sect. I3~31) . 
the normalized logarithmic conversion factor for main- 
sequence galaxies, log ( Q co/( QCO ) MS ) j can no longer be 
written as a function / ( sfr /(sfr.) M s) of normalized SFR 
(or, equivalcntly, sSFR when considering a fixed bin 
of stellar mass). This was possible, however, in the 
case of SFE and / mo i. (see Sects. 15. f .21 and 15.2. 1} 
and led to a redshift- and mass-independent recipe for 
the evolution of the population average of these quan- 
tities with sSFR/(sSFR) MS - In Fig. [Up we plo{3 
(fine dashes) the predicted variation of the median aco 
of main-sequence galaxies for three stellar mass bins 
(M*/M =5xl0 9 , 5xf0 10 and 5xfO n ) and three dif- 
ferent redshifts (.2 = 0, f, 2). Due to the higher or- 
der (s)SFR-terms in eq. [55] these trends are no longer 
redshift- and mass-independent; while aco _va hies for 
SFGs vary little across the main sequence in our high- 
est mass bin (M*/Mq = 5x 10 11 ), evolution by approx. a 
factor two is predicted between ±4 ctms for stellar masses 
M* ~ 5x 10 9 Mq. For the mass and redshifts consid- 
ered here, the mass-dependence of the average trends 
at fixed redshift is more pronounced than the redshift- 
dependence at fixed mass. Note that although we as- 
sume the relation between boost and aco-dccrement 
for starbursts (see eq. |2"o]) to be independent of red- 
shift and stellar mass, the predicted average aco-trends 
for starbursting sources nevertheless vary with stellar 
mass and redshift. This is a consequence of the mass- 
and redshift-dependency found for the 'parent' main- 
sequence population. For example, the fact that at low 
masses conversion factors are predicted to rise across 
the main sequence implies that at fixed boost-dependent 
aco-decrement the aco of high-sSFR starbursts will be 
higher than for the highest mass bins where aco _va hies 
of normal galaxies on the main-sequence are expected to 

19 All predictions shown in Fig. 114b and Figs. 1151 to [T71 assume 
the 'direct' boost-function and the corresponding best-fit value of 

7aco • 



be virtually constant. As a final comment on the de- 
scription of aco _vai "i a tions for starbursting sources we 
should point out that in practice the conversion fac- 
tor cannot decrease indefinitely (as formally implied by 
eq. [2T| but that optically thin CO line-emission sets a 
lower limit (see Fig. IT4"b). Assuming local thermal equi- 
librium and a gas temperature of 40-60 K for starburst 
sources, we estimated aco in the optically thin approxi- 
m ation using standard for mulae (see, e.g., Appendix Al 
in lBrvant fc Scovillelll996l) and obtained values ranging 
between 0.45 and 0.75 M Q (Kkm/spc 2 ) -1 . Values at the 
lower (higher) end of this range are generally predicted 
for lower (higher) redshift sources due to the evolution 
of the temperature of the cosmic microwave background, 
and with an additional contribution from the likely quite 
mild evolution of the d ust temperature in starbursts (e.g. 
iBethermin et al.ll2012ft . In relative terms, at all redshifts 
< z < 2.5 this is about 10%-20% of the aco-values ex- 
pected for massive main-sequence galaxies if their con- 
version factors also increase with redshift because of the 
general evolution of the population toward lower metal- 
licity. 

The median aco of the total population - computed 
in analogy to the bulletized procedure sketched in Sect. 
15.1.21 - is plotted with thick, solid lines in Fig. [T4"b . The 
exact shape of the transition between main-sequence and 
starburst-regime depends both on the assumed scatter 
of the FMR a nd the dispersion of aco at fixed rnetal- 
licity (see, e.g. IMannucci et al.H2010L iGenzel et aL1l2012L 
iSchruba et al.l l2012t) . For the present case we assume 
these to be 0.05 dex and 0.1 dex, respectively, which leads 
to a step-like decrease by about a factor 2-3 at an sSFR- 
excess sSFR/ (sSFR)ms ~ 3-4 with respect to the main- 
sequence average. Just as for the predicted slope of the 
normalized aco y s. sSFR relations for normal and star- 
burst galaxies, this jump changes with redshift and stel- 
lar mass. In the next section we thus provide a more 
complete mapping of expected aco _var i a tions for SFGs. 

5.3.2. aco-' predicted variations in the SFR-M* plane 

To conclude this section on empirical recipes for the 
CO-to-H2 conversion factor aco we map its predicted 
variation within the M* vs. SFR plane for three dif- 
ferent redshift bins: z = - Fig. [T5] z = 1 - Fig. [TBJ 
z = 2 - Fig. [17] We do this explicitly because the map- 
ping of metallic i ty into the M* vs. SFR plane following 
IMannucci et al.l (|2010t) is such that the variation of the 
metallicity-dependcnt aco WG adopt for 2-SFM frame- 
work is not self-similar (i.e. independent of stellar mass 
and redshift as was the case for SFE and / mo i.; see Fig. 
[T0|) and hence cannot be represented with a single, sSFR- 
dependent recipe. 

The individual panels of Figs. [13] to [T7] show the varia- 
tion of aco in M*-SFR space for the total SFG pop- 
ulation, for main-sequence galaxies and for starbursts 
(clockwise from lower left corner). We have superim- 
posed contours of constant aco an d in particular in- 
dicated the isolincs for Milky Way-like and ULIRG-likc 
conversion factors with a bold red and blue line, respec- 
tively. Due to the increasing normalization of the main 
sequence with redshift, metallicities at fixed stellar mass 
decrease with redshift (this reflects the well-established, 
measured evolution of the mass-metallicity re l ation, e.g. 
iKobulnickv fc Kewlevl 12004 lErb et al.l [20061 iLiu et al.l 
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Fig. 15. — : Predicted variations of the CO-to-H2 conver- 
sion factor, aco(f° r the J = 1 — > transition and assum- 
ing the 'direct' boost-function), in the SFR- A/* plane 
for main-sequence galaxies (panel a), for starbursting 
galaxies (panel b) and for the combined population of 
normal SFGs and starbursts (panel c). aco _v ariations 
are mapped within ±5 times the dispersion (cms) of the 
main sequence for the normal galaxies (panel a) and be- 
tween -5 <tms and arbitrarily high sSFR-excesses for star- 
burst galaxies and the total SFG population (panels b 
and c). Lines of constant aco = 0.8, 2, 3.5, 4.4, 6.5, 10 
and 16 are superimposed in blue, purple, magenta, red, 
orange, yellow and white (standard values of a ULIRG 
and Milky Way conversion factor - aco = 0.8 & 4.4 - 
are highlighted with bold lines) . Values plotted in panels 
(a) to (c) represent the median for the respective (sub- 
population. Panel (d) shows the scatter of aco (visual- 
ized here by the interquartile range IQR of logarithmic 
aco-values) around the average trends for the total SFG 
population (cf. panel c) at stellar mass M+/Mq = 5xl0 9 
(black), 5xl0 10 (light grey) and 5xlO n (dark grey). Ver- 
tical green lines correspond to the lines of constant sSFR 
plotted in panels (a) to (c). The dispersion rises strongly 
over a fairly small range of sSFR where P B ~ 50% (see 
also Figs. [Sland [Tib). 



I2008f ). While this evolution to lower enrichment is 
expected to be quite strong at the smallest stellar 
masses plotted in these figures (M* ~ 10 9 M©) the evo- 
lution is less strong for those galaxies of stellar mass 
AI± ~ 3x 10 10 M© that contribute most to the cosmic 
SFR-density over the redshift range considered here (e.g. 
iKarim et ai1l2011[ ). As a consequence, we expect that the 
conversion factor of such galaxies remains quite similar 
to the classic Milky Way-value of 4.4M (Kkm/spc 2 ) -1 
over the range < z < 2. Specifically, for a galaxy of stel- 
lar mass M* ~ 3x 10 10 M© that is located directly on the 
average main-sequence locus, the recipes developed in 
Sect. 15.3.11 predict aco — 4-1 in the local universe and 
Qco- 5.3 at z^2. Note that this prediction is not the 
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Fig. 16.—: As in Fig. [HI but for z=l. 



coincidental result of our choice of a specific slope and/or 
normalization of the relation between aco and metallic- 
ity in eq. [2] It would also hold for any of the other re- 
lations shown in Fig. [2] as these all attain quite similar, 
Milky Way-like conversion factors around solar mctallic- 
ity. 

For the variations of the conversion factor of starbursting 
sources in the M* vs. SFR plane (upper right corner of 
Figs. IT5"1 to ITT1) two different regimes exist. At high stel- 
lar masses, the aco-values of the parent, main-sequence 
population vary little (both across the main sequence 
with a given bin of stellar mass and between stellar mass 
bins); the boost-dependent aco _ decrement alone hence 
determines the value of starburst conversion factors. As a 
consequence, lines of equal starburst-aco are nearly par- 
allel to the main-sequence locus. At low stellar masses, 
lines of constant starburst-aco run nearly perpendicular 
to the iso-lines on the main-sequence locus. This is due 
to the rapid variation of aco for normal galaxies, which 
has the effect that ever higher boost-amplitudes are re- 
quired in order for starbursts of successively lower stellar 
mass to reach equal absolute values of aco (e.g. the stan- 
dard local ULIRG-value 0.8 M© (Kkm/spc 2 )" 1 ). Start- 
ing at an (s)SFR-excess of about +3ctms, starbursts be- 
gin to dominate the main-sequence population by num- 
ber. The transition between the main-sequence locus and 
the starburst-dominated part of M*-SFR space is char- 
acterized by both a sudden drop of the average aco (see 
Fig. Q3b) and an abrupt increase of the dispersion of 
aco, which is a result of the heterogeneous mixture of 
starbursting and high-sSFR main-sequence galaxies in 
this transition region. This is illustrated in the inset in 
the lower right corner of Figs. [T5l to [T7l where we plot the 
evolution of the interquartile range of aco-values mea- 
sured in the total SFG population (i.e. including both 
starbursts and normal galaxies). On the main sequence 
(sSFR/(sSFR)ms ^ 3) the scatter in aco is caused by 
the mctallicity dispersion of the of FMR at fixed M* and 
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Fig. 17.—: As in Fig. El but for z = 2. 

SFR plus the dispersion of aco at fixed metallicity and 
is hence relatively small. The aco-scatter for starburst 
galaxies is larger than that on the main- sequence locus 
because it reflects both the dispersion of aco at fixed 
M* and SFR on the main sequence, and the fact that 
the shape of the boost-function implies that starbursts 
at a given sSFR/(sSFR)ms have been boosted to higher 
(s)SFR starting from a range of positions on the main 
sequence. 

6. DISCUSSION: TOWARD A SIMPLE 
DESCRIPTION OF MOLECULAR GAS IN 
STAR-FORMING GALAXIES 

6.1. The boost-function: astrophysical context and 
limitations 

6.1.1. Star formation enhancements in simulations and 
observations 

Based on a suite of simulated interacting galaxies (with 
comparable masses and a representative range of both or- 
bital c onfig urations and mor p hologi es) IDi Matteo et al.l 
(pOTJ7t) and IDi Matteo etaLl ((20081 ) derived the SFR- 
cvolution of major mergers as compared to the evo- 
lution of identical, isolated galaxies. In Fig. @] we 
plotted the ma xi mal SFR-enhancements reported in 
IDi Matteo et al.l (|2008D . averaged between the Trce- 
SPH (smoothed particle hydrodynamics) simulations and 
grid-based N-body simulations carried out by these au- 
thors. Both the 2-SFM boost-function and the distri- 
bution of simulated SFR-enhancements have a clearly 
defined peak. The exact position of this peak (which 
corresponds to (x)bk m the 2-SFM formalism) in the 
simulations depends on the gas content of the galaxies. 
In 'local' disk galaxies (simulated total gas fractions be- 
tween 10 and 30%) the most frequently encountered max- 
imal SFR-excess is approx. a factor of three, while for 
gas-rich simulated galaxies (/ mo i. ~ 50%) reminiscent of 
high-z disks it is twice as large, mainly because these tend 
to become Jeans-unstable and form dense gas clumps 



when perturbed. The 2-SFM boost-function peaks at a 
roughly four-fold SFR-enhancemcnt but it should not be 
directly compared to the simulation results because not 
all starbursts that contributed to the shape of the un- 
der lying (s)SFR-distribution in iRodighiero et all (|2011[ ) 
can have been 'caught' at the peak of the starburst ac- 
tivity. Wit h respect to the results of IDi Matteo et al.l 
(|2007ll2008D . observable distributions of SFR-boosts for 
interacting galaxies will likely be modified if minor merg- 
ers and fiy-bys, as well as the relative timing of star- 
burst events, are accounted for. If SFR-enhancements in 
merger-driven starburs ts depend on the mass ratio of the 
galaxies involved (e.g. ICox et aTl 120081 ). then the boost 
distribution including minor mergers should be broader. 
Likewise, fiy-bys and asynchronous burst-activity plausi- 
bly shift and skew the boost-distribution to lower SFR- 
enhancements . 

Although the distribution of SFR-enhancements caused 
by interactions between galaxies is poorly constrained, 
observationally, there have been numerous studies aimed 
at quantifying the integrated contribution of excess star 
formation associated with mergers /inter actions to the 
cosm i c SFR density at redshi f ts z < 2 (e.g.lRobaina et al.l 
[20091 IKampczvk et^[20TTl lKaviraietal.1 120131 ). One 
attempt is the recent determination of the distribu- 
tion of SFR-enhanc ements in SDSS galaxy pairs by 
iScudder et al.l (|2012f ) which we plot in blue in Fig. 0] 
Compared to the b oost distribution i n the simulations 
of major mergers in IDi Matteo etH] (|2008l) it is indeed 
displaced to systematically lower SFR-enhancements. In 
addition to the expected shifting and skewing, the fact 
that the SDSS pair sample does not include mergers in 
which final coalescence has already taken place implies 
that their distribution of SFR-enhancements represents 
a lower observational limit to the total local SFR-excess 
distribution caused by inter actions. The simulation s of 
IDi Matteo et al.l (|2007D and IDi Matteo et al.l (|2008f ) en- 
compass a broad variety of orbital configurations of merg- 
ing galaxy pairs and have been statistically weighted to 
reflect the dependence of the collision rate on the rel- 
ative velocities and impact parameters. However, they 
do not provide information on the contribution of minor 
interactions, nor are they carried out in a fully cosmo- 
logical framework that accounts for, e.g., the preferen- 
tial alignment of ga laxies in different locations within 
the cosmic web (e.g. lHahn et al.ll2010D. This additiona l 
step was taken in simulations by iHopkins et al.l (|2010D , 
such that their spectrum of merger-induce d surplus SFR 
should be similar to the measur ement of IScudder et al.l 
(|2012t) . IHopkins etaLl (120101) incorporated the re- 
sults of their own h igh-resolution merger simulations 
(jHopkins et al.l [2009T) in a cosmological (dark-matter) 
framework and were able to predict SFR-distributions 
of s ecularly evolving and starbursting galaxies (see Fig. 
7 in IHopkins et al.ll2010l ) that, in qualitative terms, re- 
semble the split into main-sequence and starburst ac- 
tivity we proposed in S12 and which seem broadly con- 
sistent with th e expected modificatio ns to the boost- 
distributions of IDi Matteo et all (|2008l ) discussed at the 
end of the last paragraph. For galaxies with a stellar 
mass of 10 11 Mq, where we can compare with our own 
double log-n ormal decompo s ition a ccording to eq. [SJ the 
approach of IHopkins et al.l (120101 ) predicts (1) a SFR- 
boost distribution for starbursts that is broader and more 
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skewed to low boosts, and (2) typical SFR-enhancements 
that are smaller (for simulated z ~ and 2 mergers) in 
comparison to both the outcome of t he major merger 
simulations of iDi Matteo et all (|2008fl and also the 2- 
SFM boost-function. 

To infer that galaxy-galaxy interactions cannot be the 
sole trigger of starburst activity, based only on the mis- 
match between the 2-SFM and measured or simulated 
boost distributions, would, however, be premature. The 
approach of decomposing an sSFR-distribution into two 
components - as done in Sect. 14.21 - leads to inher- 
ently poor constraints on the shape of the boost-function 
at low boosts (see cross-hatched area in Fig. 2]) since 
galaxies with small sSFR-enhancements blend in entirely 
with the main-sequence population. This 'maximiza- 
tion' of the main-sequence contribution in eq. [5] will 
hence cause a truncation of the lower part of the boost- 
function. On the other hand, the distinction between 
secularly-evolving and only weakly starbursting systems 
itself is not clear-cut because minor merger events occur 
frequently. If all galaxies that have experienced mini- 
mal boosting are regarded as starbursts, then a boost- 
distribution that goes to zero at a (s)SFR-excess of zero 
(boost = 1) is unrealistic. The 2-SFM boost- function 
should thus best be viewed as the signature of strong 
boosting where a significant fraction of the ISM fuels 
starburst activity. Despite this limitation it is interest- 
ing that the peak position of the 2-SFM boost-function 
at an excess (s)SFR of a factor of four corresponds ex- 
actly to the a verag e SFR-boost that was measured by 
iHwang et al.l (|2011l) for FIR-selected galaxies at z = 
and z = 1 undergoing an interaction with a late-type 
neighbor. In an alternative measurement of interaction- 
induced SFR - boost ing in local late-type galaxy pairs, 
IPark fc Choil (f2009h found an increase of the equivalent 
width of the Ha-line by an identical factor four when 
the two galaxies were separated by less than 1% of the 
virial radius of the companion's halo. This constancy of 
the average SFR-cxcess is reminiscent of the evidence for 
only mild evolution of the shape of the boost-function 
we presented in S12. Furthermore, it is worth noting 
that notwithstanding the uncertainties concerning the 
shape of the lower en d of the 2-SFM boost-function, 
iBethermin et al.l (|2012t ) successfully used it as the basis 
for matching observed IR source counts. Their analysis 
was a good test of the viability of the 2-SFM descrip- 
tion of starbursts since it employed different IR SEDs 
for main-sequence galaxies and starbursts. 

6.1.2. Link to dark matter 

Merging and star-formation activity - even in the "sec- 
ular" mode - reflect the accretion of dark matter (DM) 
and the primordial gas bound to the DM haloes. It is 
thus interesting to check whether there are clear simi- 
larities between the distribution of SFRs of galaxies of 
a given mass and the a ccretion of DM ont o the corre- 
sponding parent haloes. iDekel et al.l (|2009f ) determined 
the DM infall rates at the virial radius of >100 simu- 
lated DM haloes with mass 10 12 M at z~2.5 which 
typically host ~10 n M® galaxies. For these systems 
(w hich have masses comp arable to the galaxies used 
by iRodighiero et all 120111 to construct distributions of 
sSFR) the DM accretion spectrum shows an extended tail 
of high accretion rates which is dominated by "major" 



merging activity where the mass ratio between accreted 
and parent DM halo is fairly high. It is obviously tempt- 
ing to associate this feature of the DM ac cretion rate 
distrib utions in the simulations presented in IDekel et al.l 
(|2009t) to the tail of excess SFRs contributed by star- 
burst galaxies while the smooth accretion would then 
fuel the sustained secular mode of star formation that is 
characteristic of SF Gs on the main seq uence. This was 
already proposed by IDekel et al.l (|2009f ). who also point 
out that in this context the abbreviation "SFG" could 
legitimately stand for "stream-fed galaxy". In analogy 
to our split of the sSFR-distribution in eq. [51 T. Goerdt 
et al. (in prep.) have decomposed the DM accretion rate 
distribution for such z = 2.5 DM haloes with virial mass 
10 12 Mq into two log-normal contributions and find that 
the one shifted to high accretion rates contributes ap- 
prox. 10% to the infalling mass budget. This value is 
strikingly similar to the 14.2l^g% (68% confidence lim- 
its) we inferred in S12 for the contribution of burst-like 
star formation to the total SFRD at z = 2. 

6.2. Universal star formation laws and the distribution 
of galaxies in the Schmidt- Kennicutt plane 

The tightness of the star formation law (dispersion 
~0.2 dex or less than a factor of two) which we found us- 
ing our newly 'homogenized' literature data in Sect. 13.21 
is remarkable and points to a very direct and apparently 
ubiquitous link between the global molecular content of 
galaxies and how much of it is being converted into stars. 
It is akin to stating that, for normal disk galaxies out to 
at least z ~ 2.5, once the SFR has been measured the size 
of the associated molecular gas reservoir can be inferred 
with high accuracy, and vice-versa. In our calibration of 
the integrated S-K law in Sect. [3J we adopted a statis- 
tical (metallicity-dependent) estimate of the CO-to-Eb 
conversion factor, qcOi f° r 90% of the normal galaxies 
in our reference sample when we were translating the 
observed correlation between Lib. and L' co to a more 
physical relation between SFR and M mo i.. While using 
an average «co and neglecting the associated scatter in 
principle artificially reduces the dispersion of the S-K 
law, it is not inconceivable that a dispersion in aco has 
produced the observed width of the Lyr, vs. L' co re- 
lation and that the underlying star formation law might 
be intrinsically even tighter than the 0.2 dex we measured 
here. To truly test the universality of the 'normal'-galaxy 
S-K law there are at least two complementary ways for- 
ward. On the one hand, it will be important to compile 
samples of galaxies for which star formation and gas esti- 
mates rely on strictly identical tracers (e.g. the ground- 
state transition of 12 CO). On the other hand, it may 
prove worthwhile to assess in detail (e.g. by means of 
a full sampling of the SLEDs of molecular gas tracers) 
how star-formation proceeds in different phases of the 
ISM. Having the capability of doing this in a resolved 
fashion, we also have the potential to reveal what causes 
the mild SFR-dcpendcncc of SFE, which manifests it- 
self as a non-lin ear slope of the integra ted S-K relation 
(SFRocM^ ). ISaintonge et al.l (|2012fl have proposed 
that the rise of SFE across the main sequence is due to 
morpho logical 'stabilization ' of the ISM in bulged galax- 
ies (see iMartig et al.l I2009D , which are more abundant 
on the lower part of the main-sequence locus. Resolved 
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studies of the star formation law will also be able to re- 
veal whether a similar mechanism is responsible for the 
SFE-increase with rcdshift in massive galaxies; while it 
could be due to the inc reasing absence of bulged galax- 
ies at high rcdshift (e.g. lOesch et al.ll2010h . it seems just 
as plausible that star formation in an increasingly tur- 
bulent medium including massive star-forming clumps 
would proceed in a more efficient fashion. 
The variations of SFE with SFR among normal galax- 
ies are small compared to the strong SFE-enhancemcnts 
that are observed in starbursting systems. In Sect. 15.1.21 
we introduced an empirical, supra-linear scaling between 
the SFE- and (s)SFR-increase during starburst episodes. 
This relation can be understood very intuitively by the 
balance between the three quantities involved: SFR and 
M*, which grow in the burst-phase and M mo i, which de- 
creases as gas is converted into stars. SFE, as the ra- 
tio between SFR and M mo i., thus inevitably increases 
more strongly than sSFR. An immediate consequence 
of this is that we expect a more spread-out distribu- 
tion of star formation efficiencies than is observed. In 
Fig. [9] we illustrate how, in the 2-SFM framework, the 
overlapping sSFR-distributions of normal galaxies and 
starbursts move apart into a more clearly double-peaked 
SFE-distribution. Discrete recipes for the assignment of 
CO-to-H2 conversion factors to normal galaxies and star- 
bursts arc thus not the only way to obtain a bimodal 
distribution of galaxies in the S-K plane; this can also 
be achieved with a more physical, continuous descrip- 
tion of star-formation in starbursts. Width and depth 
of the trough we predict between the 'sequence of disks' 
and 'sequence of starbursts' in the S-K plane depend on 
the shape of the l owe r end of the boost-function (see 
discussion in Sect. 16. ip . Our model of an unbiased pro- 
file through the S-K plane at fixed gas mass does high- 
light, however, that observations of large cosmological 
volumes are necessary to fully sample the actual distri- 
bution of galaxies with respect to SFE: the relative am- 
plitude of the SFE-distributions of starbursts and main- 
sequence galaxies in Fig. [9] is expected to be a factor of 
30 and the contrast between the peak of the starburst 
distribution and the trough merely a factor two. A first 
attempt to construct a representative sampling of S-K 
space using the COLD -GASS survey was presented in 
iSaintonge et al.l (|2012t see their Fig. 6b) and demon- 
strated just how insignificant starbursts are in determin- 
ing the shape of the star formation law for the bulk of 
the population. 

6.3. The consumption of gas reservoirs during 
starbursts 

As a consequence of the high efficiency with which gas 
is converted to stars in starburst episodes, the gas reser- 
voir in the host galaxy is used up more quickly than 
it can be replenished by accretion from the intergalac- 
tic medium. In Sect. 15.2.11 (see Figs. [TU] and [11]) we 
showed that the overall gas fractions (i.e. gas fractions 
taking into account the molecular and stellar mass con- 
tent throughout the whole starbursting galaxy) of star- 
bursts in our reference sample are indeed in general lower 
than the average gas fraction of galaxies which reside on 
the main sequence. We can use the 2-SFM description 
of starbursts to explicitly calculate how we expect the 
gas content of starbursting galaxies to change once it 



has left its main-sequence state. Given that each of the 
starbursts in our calibration sample was observed in a 
different stage, this should be viewed as a comparison be- 
tween the gas fraction prior to the onset of the burst and 
the gas fraction which would be measured approximately 
half-way through the starburst event. If we consider the 
gas-to-stellar mass ratio /i mo i. = M mo i./M*, the relation 

between pre-burst (^™f bmst ) and mid-burst gas content 
takes a simple and only boost-dependent form (sec Fig. 
H8b): 

/W./M£o e iT bUr8t = (boosWR) 1 -^ . (29) 

Here we have used that /i mo i. = sSFR/SFE 
and that the sSFR and SFE of the star- 
burst are sSFR /sSFR prc ~ burst = boost s sFR and 
SFE/ SFE prc-bu ra t _ (boost s sFR.) 7SFE , respectively. In 
Fig. [TSh we plot the variation of the typical ratio 
between mid-burst and pre-burst gas-to-stellar mass 
ratio as a function of sSFR-excess. This average trend 
is the result of pairing up each point on the dotted 
curve for the evolution of /U mo i. (sSFR) for starbursts in 
Fig. HUb with a position on the corresponding relation 
for normal galaxies (dashed line in the same figure) 
by means of the boost-value at the peak of the 

sSFR-dependent boost distribution (see Fig. [5] and 
Sect. 14.2. 3|) . Note that this calculation assumes the 
stellar mass in the starbursting galaxy and its pre-burst, 
main-sequence state to be equal. By neglecting the fact 
that stellar mass has been added to the system during 
the first phase of the burst, our estimate of the ratio of 
mid-burst to pre-burst gas fraction effectively represents 
an upper limit. This simplification also makes the 
average trends in Fig. [18] independent of stellar mass 
and redshift, because all dependence of the absolute 
value of the pre-burst gas fraction on these two factors 
(see e.g. Fig. [13]) is eliminated. The hatched/shaded 

regions straddling the median trend for /U mo i.//^oi bmst 
reflect the 1 a uncertainty on the relation between 
average boost and sSFR (see also Fig. [S^). 
Based on the theoretical understanding derived from the 
2-SFM approach, we are for the first time also able to 
infer - in a statistical sense - the pre-burst gas fractions 
of the starburst galaxies in our reference sample, i. e. of 
eight local ULIRGs from lDownes fc Solomon| (Il998f) and 
of three high-redshift starbursts studied bv lMagdis et all 
(|2012bl) . We do so under the same assumptions as 
already used to derive the theoretical curve discussed 
above and superimpose our estimates on the 2-SFM 
prediction in both panels of Fig. [18] Given that the 
typical sSFR-excess of our reference starbursts is about 

a factor ten, their median /^moi./(A i mor bmSt ) °f ~0.35 
is in quite good agreement with the 2-SFM prediction 
for the scenario of the 'direct' boost-function (black 
line in Fig. [TEk). For the most common starbursts, 
which have a (s)SFR-boost equal to four, the molecular 
mass-to-stellar mass ratio half-way through the burst is 
expected to lie between 50% (30%) of its initial value 
for the direct (merger-corrected) boost-function. 
It is interesting to explicitly compare our constraints 
on the gas fraction decrease during the starburst phase 
with that expected in the case that the starburst is 
triggered by a major merger. We approximate the 
SFR-evolution during the interaction-induced burst 
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Fig. 18. — : (a) Predicted ratio between the molecular gas-to-stellar mass ratio of starbursts, /i mo i. = Af mo i./Af*, and 

the average expected mass ratio (Mmoi~ UrS *) P r i° r to the onset of the burst (i.e. during a phase of secular growth on the 
main sequence). The variation with sSFR reflects the different sSFR-dependence, for normal galaxies and starbursts, of 
molecular gas-to-stellar mass ratios predicted by the 2-SFM description (see Fig. [TUb ; dashed and dotted lines, resp.). 
All sy mbols are as in Fig . El Th e pre-burst gas-to-stellar mass ratios of starburst galaxies from iDownes fc Solomon! 
(|1998[) and iMagdis et al.l (|2012bf ) were inferred assuming that these sources experienced the median (s)SFR-boost 
expected for sources with the same sSFR-excess (see Fig. [5]). 

(b) Dependence of the ratio /imol./(Mmoi~ UrS *) on the sSFR-boost of starburst galaxies (see eq. [29]) . The black (grey) 
line shows the 2-SFM prediction for the average trend for the direct (merger-corrected) boost-function. 



by a top-hat function such that, at a time i™ c 2 rscl 
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and rearrange terms to obtain an expression for the ratio 
of the pre- and mid-burst gas fractions: 
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Here we used that the total mass (Af mo i + Af*) stays 
constant and that the SFE before and during the star- 
burst are related by the power-law in cq. [23l With 
a typical boost of approx. a factor six and 7sfe = 1-5 
(direct boost-function) we find /moi.// m o C i~ burst ~ 0.6 for 
r prc-burst ^ -y Qy r g ag d e pi e yon time scale of main- 
sequence galaxies) and ^ e 2 rser « 50 Myr (we take this 
number to be about half the time for which interaction- 



induced star formation is sustai ned in numerical sim- 
ulations of galaxy mer gers, e.g. iDi Matteo etal] l200l 
iBournaud et al.l [201 lbf) . Accounting for the mass- and 
rcdshift-dependence of the conversion^ between relative 
gas fractions and gas-to-stellar mass ratios, this corre- 
sponds to values of Mmoi./MmoiT UMt m the range of 0.4 
to 0.5. Obviously, the simple calculation leading up to 
cq. [30] will in reality be complicated by, e.g., gas los s 
and heating in merging systems (e.g. ICox et alll2004|) . 
the modified balance between the atomic and molecular 
hydrogen phase in dense, turbulent media, and IMF vari- 
ations as have been proposed for starburst regions (e.g. 
Baugh et a l. 2005, Papadoooulos et al~ll201lL but see also 



Tacconi et al.ll2008l iHavward et al.ll2013[ ). Taken at face 



value, the reasonable consistency of the estimates of the 
gas fraction decrease as per eqs. |2"91 and |3"01 may indicate 
that neither of these three factors plays a major role (or 
that these competing effects compensate each other) . 
Systematic comparisons between the molecular gas frac- 
tions of normal and starbursting galaxies will reveal 

20 The relative gas fraction and gas-to-stellar mass ratio of the 
pre- and mid-burst state are related by 

/ . pre — burst 

r I /,pre — burst / pre — burst / ^mol. 

Jmol./j mo l. — Mmol. /Mmol. I , T 

\ J- + Mmol. 



For the massive SFGs discussed here, an initial gas-to-stellar mass 
ratio urst of ~5-10% and 50-100% is expected at low and 

high-redshift, respectively. The term in brackets should thus vary 
between roughly 1.1 and <2. 
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whether the trends we proposed based on our small sam- 
ple are robust. Further tests of the 2-SFM framework 
will now be discussed in Sect. 16.41 



6.4. Observational validation of assumptions and 
predictions made by the 2-SFM approach 

The 2-SFM framework as we have developed it so far 
has produced a remarkably simple description of SFGs 
over the last 10 Gyr. One may legitimately wonder 
whether this simplicity is the true imprint of fundamental 
laws that govern galaxy formation in a cold dark mat- 
ter Universe or the outcome of an incomplete or selective 
view of the star-forming population due to observational 
limitations. The answer to this question depends to some 
extent also on the scope of any investigation. The occur- 
rence of extreme behaviour in rare outliers or small scale 
processes with little impact on global system properties 
- while relevant for a complete understanding of all com- 
plex aspects regulating star formation - does not imply 
a general inadequacy of a simpler approach, as we have 
been advocating here, which aims to provide a panoramic 
treatment. Further confirmation of the validity of the 2- 
SFM description will instead involve both (a) revisiting 
some of its key ingredients and (b) testing its predictions. 
Concerning point (a), the main focus should lie on veri- 
fying our hypothesis that galaxies with stellar mass sig- 
nificantly below M±/M Q = 10 10 follow the same relations 
that were calibrated on galaxies which are more massive 
than this threshold. The universality of the S-K law, for 
example, will soon be routinely tested with ALMA down 
to low stellar masses and out to high redshift by targeted 
observations of lensed galaxies. With deeper follow-up 
of molecular transitions, it will also be po ssible to iden- 
tify e volution in the normalization (e.g. iTacconi et al.l 
120131 ) and curvature of the star formation laws. A sec- 
ond assumption of the 2-SFM framework is that the dou- 
ble log-normal decomposition of the sSFR-distribution 
is applicable also at z ^ 2. To ascertain this, trac- 
ers of dust-obscured star-formation are indispensable as 
extinction-corrected SFR-mcasurem ents underestimate 
the true SFR of dusty starbursts (e.g.lHughes et al .111 9981 
Trentham et al.1 119991 iBuat et al.1 120051 iChapman et al.1 

i]F~ 



20051 iDaddi et al.l l2007al iCasev et aTTl2013[ ) and place 
these on the locus of the star-forming main sequence. 
Obtaining good statistics on the rare starbursting sources 
(comoving number densities are of order 10 -5 Mpc~ 3 ) 
at the high-end tail of the sSFR-distribution hence re- 
quires a combination of wide-area IR or radio surveys 
with complementary, deep optical or UV data. A non- 
universality of the double log-normal decomposition in 
cq. [5] would introduce more variation in the simple 
SFE vs. sSFR-excess relations, etc. than is currently 
suggested by the fairly limited data. Any evolution in 
the sSFR-decomposition into normal galaxies and star- 
bursts would imply a more complex, redshift-dependent 
behaviour of average scaling relations in the space of nor- 
malized molecular gas properties than is shown in Fig. 
ITOl A more fundam ental question is whethe r the z ~ 2 
sSFR-distribution of [Rodighicro et al.l (|2011l ). on which 
we perform the decomposition to begin with, is accu- 
rate. Little is known about the lower tail of the sSFR- 
distribution, but it is unlikely that low-sSFR outliers to 
the main sequence should be responsible for a signifi- 
cant amount of star formation activity (see S12 and ref- 



erences therein). In the absence of a single star formation 
tracer to map out the distribution of ga laxies in the SFR- 
pl ane, the two-pronged approach of lRodighiero et all 
(| 2 1 If) for reconstructing it with two different diagnos- 
tics (IR emission for dust-obscured galaxies and UV- 
emission for the bulk of the main-sequence population) 
relies on the consistency of the associated star formation 
estimates. Extinction-corrected UV-fluxes and IR mea- 
sure ments at z ~ 2 are kn own to agree in an average sense 
(e.g. IDaddi et al.ll2007al ). but the dispersion about the 
mean extinction-correction could potentially contribute 
to the observed scatter of the main sequence of SFGs. 
The analysi s of th e main sequence at 0.5<z<1.3 by 
I Salmi et al.l (|2012[) suggests that at least at these red- 
shifts th e dispersion of the sequ ence is mainly intrinsic, 
bmce m iRodighiero et al.1 (|2Qlll ) the sSFR-distributions 
of galaxies with M* > 10 11 Mq are identical when com- 
puted with UV- or IR-cmission, this seems to hold for 
the high-mass end of the z ~ 2 main sequence as well. 
Concerning point (b), the analysis of this paper produced 
predictions that will be tested in future CO follow-up 
observations. With a good sampling of the transition 
region between main sequence and starburst galaxies in 
the SFR-M* plane (sSfr/ (s 

SFR) MS € [3, 5]), these observa- 
tions will quantify the scatter of, e.g. SFE and determine 
whether it is indeed larger than elsewhere, as is expected 
for a heterogeneous mixture of normal galaxies and star- 
bursts. While we predict such an increased dispersion to 
be measurable even using direct observables, e.g. L' co 
and I/ir, an estimate of the CO-to-E^ conversion fac- 
tor aco is necessary to calculate actual values of SFE 
and f gas . Our predictions for the variation of aco m the 
SFR-M* plane in Sect. 15.3.21 are in principle testable, 
but obtaining high-confidence measurements of aco will 
remain a challenging task that is best tackled using dif- 
ferent, complementary strategies in parallel. The gas - 
to-dust ratio t echnique empl o yed bvlLerov et al.l (l20Tll) , 
IMaedis et all (IMTlf I2012bfl . IMaenelli et al.l (|2012h is 
powerful, in that it can provide constraints on the con- 
version factor for large data sets. However, the large 
scatter in measured gas-to-dust rat ios in local, low- 
metallicity and low- mass galaxies fe.g.lDraine "&I1I2007L 
IGalliano et all 120081 IGalametz et all 120 lit Remv et al, 
in prep.) indicates that this method becomes highly in- 
accurate for z > 3 galaxies and galaxies with stellar mass 
M* <C 10 10 M Q . An alternative approach is to interpret 
dynamical constraints from CO line profiles in the con- 
text of numerical sim ulations to infer aco as proposed by 
IDaddi et all (|2010al ). However, even barring the system- 
atic uncertainties on model dark matter distributions, 
the application of this method to large data sets may be 
impracticable as it requires high signal-to-noisc data and 
a fine spectral sampling of the emission feature. 
Finally, we note that the (s)SFR-boost of starbursts, al- 
though much less easily determined than their (s)SFR- 
excess with respect to the average of the main-sequence 
population, is in principle measurable using high-fidelity 
and ideally also spatially resolved spectroscopy. When 
compared to the output of stellar evolution models, this 
kind of data would allow a detailed reconstruction of 
the SFH of boosted sources prior to the onset of burst 
activity. It would hence also reveal whether galaxies 
that show a strong sSFR-excess are truly experiencing 
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short-term boosting of their activity at all rcdshifts or 
whether they are merely a high-intensity tail of the 'nor- 
mal' population. The episodic and merger-related na- 
ture of ULIRGs at low redshift is well-accepted (e.g. 
lSanders"fe Mirabel 119961 : and references therein) but is 
hard er to prove for starbursts in the distant Universe 
(e.g. iTacconi et al.ll2006l iDaddi et aLll2009l Hvison et al.l 
12013ft . The supporting evidence which has been accu- 
mulating in recent years, however (e.g. IR diagnostics, 
ISM temperatures, host galaxy structure and kinemat- 
ics; see also our overview in the introduction and the 
discussion in Sect. 16 -If) , is at the basis of our proposed 
split into starburst and normal galaxy populations and 
the assumption that it provides a valid description of 
the star-forming population of much of the history of the 
Universe. 

7. SUMMARY 

The 2 Star Formation Mode ('2-SFM') framework pro- 
vides a conceptually simple and self-consistent scheme 
for the prediction of basic properties of the star-forming 
galaxy (SFG) population. It relies on basic observ- 
ables - e.g. the evolution of specific star formation rate 
(sSFR) in main-sequence galaxies or their stellar mass 
(Mj.) distribution - and their mathematical description - 
e.g. the Schechter function parametrization of the stellar 
mass function or slope and normalization of the Schmidt- 
Kcnnicutt (S-K) law - to produce an analytico-empirical 
description of the statistical properties of SFGs which 
can be both predictive and help (re-)interpret existing 
measurements. A central ingredient of the 2-SFM frame- 
work is the distinction between 'normal' SFGs that re- 
side on the star-forming main sequence and starbursts 
that are are much rarer and regarded here as a 'pertur- 
bation' of the main-sequence state (see Sect. 14. ip that 
is probably dynamically induced or induced by interac- 
tions. We recently applied this approach successfully 
for the prediction of IR luminosity functions at z < 2.5 
([Sargent et al.l I2012T ) and of galaxy nu mber counts be- 
tween 24 and 1100 /jm and at 1.4 GHz (Bcth ermin et al.1 
12012ft . In this article we have investigated the observa- 
tional evidence that the molecular gas properties of mas- 
sive (M*>10 10 Mq) SFGs are amenable to a similarly 
simplified description as their IR-emission. 
We use a sample of approx. 90 normal SFGs (see Sects. 
12.1.11 and 12 . 1 .2(1 to derive scaling relations that describe 
the molecular gas content of secularly evolving SFGs at 
z < 3. When all involved quantities are normalized to 
the value a given observable takes for an average main- 
sequence galaxy, these scaling relations become strikingly 
simple (and in general also independent of redshift). In 
particular, we find that: 

1. All literature measurements of SFR and M mo \, in 
massive (M* > 10 10 M Q ) main-sequence galaxies at 
z < 3 arc compatible with the existence of a uni- 
versal (i.e. redshift-invariant) star formation law 
for such systems. This integrated S-K relation is 
slightly supra-linear (SFR cx ) and tight (dis- 
persion ^0.2 dex; see Fig. [T]). 

2. Star formation efficiency (SFE) varies very little 
across the main sequence (see Fig. [7|) while the 
molecular gas mass fractions, M mo i./M*, increase 



almost linearly with (s)SFR for main-sequence 
galaxies of a fixed stellar mass. 

3. Changes in the sSFR of main-sequence galaxies 
are strongly correlated with changes of the molecu- 
lar gas fraction, implying that both the dispersion 
of the main-sequence and the cosmic evolution of 
sSFR in general reflect variations of the gas content 
of normal galaxies (see Fig. fT2|) . 

Based on this characterization of gas in the main- 
sequence population, we are then able to predict the 
molecular gas properties of starburst galaxies which - in 
the 2-SFM approach - start out as normal galaxies that 
subsequently experience boosting to higher (s)SFRs. By 
considering the excess SFR and excess SFE of observed 
starburst galaxies with measured CO-to-H2 conversion 
factors aco, we infer that SFE grows more strongly in 
the burst-phase than SFR (see Sect. I5.1.2j) . Taking into 
account the changing, sSFR-dependent mixture of star- 
burst and normal galaxies that constitutes the total star- 
forming population, this leads to the following expecta- 
tions: 

4. Normal SFGs and starbursts are separated more 
strongly in the S-K plane than in the space of M* 
and SFR (see Fig. [9]). However, a separation that is 
as discrete as currently suggested by observations 
is not expected and is likely the outcome of the 
incomplete sampling of the S-K plane in surveys 
explicitly targeting strong starbursts and average 
main-sequence galaxies. 

5. Even if starbursts are treated as a continuous ex- 
tension of normal galaxies, with depiction times 
that decrease in proportion to their burst-related 
(s)SFR-enhancement, a nearly step-like, roughly 
ten-fold increase of the SFE is predicted at the 
sSFR where starbursting sources begin to outnum- 
ber the main-sequence population (see Fig. fTOk). 

6. A similar, albeit less pronounced step-like be- 
haviour is predicted for molecular gas fractions (see 
Fig. [10p): while these continuously rise across 
the main sequence (see point 2 above), the higher 
SFE of starbursts causes their gas fractions to de- 
crease to a value that is smaller than the average 
obser ved for a typical main-sequence galaxy. In 
Sect. I6.3l we provide recipes for how much gas frac- 
tions are expected to drop, depending on the in- 
tensity of the starburst, if - as is expected for, e.g., 
merger-induced starbursts - the time-scale for the 
exhaustion of the molecular fuel reservoir is much 
shorter than the time-scale for accretion of pristine 
gas from the cosmic web. 

Based on the systematic difference between the L' co /Lir 
and Mh 2 /£ir ratios of starbursts, we derive an empirical 
recipe for the CO-to-H2 conversion factor, acoj of star- 
burst galaxies (see Sect. 15.3. 1[) . In combination with 
an assumed metallicity-dependence of aco for main- 
sequence galaxies, we are able to predict aco-variations 
for SF Gs throughout the M* vs. SFR plane (see Sect. 
I5.3.2[) . Due to the flatness of the mass-metallicity rela- 
tion at high stellar masses, the conversion factor of Milky 
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Fig. 19.—: Redshift-dependence of the sSFR of SFGs with stellar mass M*/M Q «5x10 9 (top) and 5xl0 10 (bottom), 
as published in the recent literature (see legend; where necessary, literature values from adjacent mass bins were 
used to interpolate to the mass scales displayed here). Measurements derived based on image-stacking are indicated 
with open symbols and error bars denote the uncertainty on the sSFR-average rather than the sSFR-scatter in the 
population. Solid/dashed black lines - the best- fit evolution of the sSFR - parametrized as in eq. IA1I (see inset panels 
on lower right for the covariance between the free parameters of the fit) - and associated 2cr-errors.; light grey lines 
- sSFR-evolution in the other of the two stellar mass bins depicted in the figure, for comparison; dot-dashed line - 
evolution according to (1+z) 28 as used in S12 for the range z < 2. 



Way-mass galaxies is expected to resemble the canonical 
Milky Way- value even at the cosmic epoch when the star- 
formation history of the Universe peaked. 
The simple, analytico-cmpirical description of molecu- 



lar gas in star-forming galaxies developed in the present 
work will be used to infer the evolution of molecular gas 
mass functions and CO luminosity functions, as well as 
CO source counts in two forthcoming papers. 



APPENDIX 

LITERATURE MEASUREMENTS OF SPECIFIC STAR FORMATION RATE IN MAIN-SEQUENCE 

GALAXIES 

The locus of the main sequence is known to depend on sample selection (e.g. iKarim et al.ll20ill ) and, in particular, 
on how actively star-forming the sample under consideration is. With the aim of deriving a representative, average 
evolution we gathered measurements from several recent studies of the distribution of SFGs in the (M*, SFR)-plane that 
employed different selection criteria (e.g. different color cuts or selection by morphology, by near-IR flux/mass or by 
SFR), chose different star- formation tracers (e.g. UV, IR or radio emission) and/or adopted a variety of measurement 
techniques (e.g. individual detections vs. source stacking). By considering two separate mass scales at 5xl0 9 and 
5xlO lo M0 we obtain a constraint on the typical exponent v of the M*-dependence of sSFR, sSFRocM*". We find 
that an exponent ^~-0.2 reproduces the systematic shift between the sSFR-evolution of galaxies in the two M*-bins 
(see Fig. [IT? )). This slope agrees clo sely with the value ^ = -0.21±0.04 we adopted (based on the z^2 main sequence 
presented in iRodighiero et al.l 120111 ) in our previous publications investigating the viability of the 2-SFM framework 
(see S12, IBethermin et al.l I2012T ) . Our literature compilation covers the redshift range z < 7, with a majority of the 
measurements tracing the steep rise of sSFR in main-sequence galaxies out to z~3. At z > 4, drop-out samples 
constrain the sSFR-evolution at M* ~ 5x 10 9 M Q , but they do not contain enough high-mass galaxies to probe the 
evolution in our second, more massive bin. 

We parametrize the evolution of sSFR with a smoothly varying function of redshift with five free parameters which 
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we fit to the data in Fig. [121 

sSFR(Af„ z) = N(AQ exp ( 1+ ^ c ) , (Al) 

where 

N(M+) = N(5x 10 10 M ) 10" lo sK/[5x io 10 m o] ) 
A = 2.05±g;|g 
B = 0.16l°;J 5 7 
C= 1.54 ±0.32 

Here ^ is the slope of the log(sSFR) vs. log (A/*) relation as above and the normalization at a stellar mass of 5x 10 10 M Q 
is N(5x 10 10 M©) = 0.095^q qq3 Gyr -1 . The quoted uncertainties are 68% confidence limits as determined by a Monte 
Carlo Markov Chain (MCMC; 10 6 realizations). 

In Fig. [TO] the 1 cr-errors on the a verage sSFR-evolution a ccording to eg. lAll are m arked by dashed lines. Due to the 
abun dant data (e.g..lNoeske et al.ll2007l lElbaz et al 1 120071 iDaddi et al.ll2007bl 120091 iPannella et"aIll2009L iKarim et all 
l20Tll iRodighiero et all 120111 iWhitaker et al.l 12012ft and their generally high fidelity (error bars on the individual 
literature measurements span the statistical uncertainty on the mean rather than the population dispersion), formal 
uncertainties at z < 2 are small but increase steadily thereafter, reflecting the much sparser data at the highest redshifts. 
The growing formal errors do not include the systematic evolutionary uncertainties at z > 3, where the sSFRs of drop- 
out ga]BxiBshaye_been su bject to frequen t revis ion on an almost yearly bas is. Initial mea s ureme nts at 4 < z < 8 
by iStark et ahl (j2009l ) and [Gonzale z et al.l (|2010T ) - subsequently modified by iBouwens et al.l (|2012l ) to account for 
dust -extinction - suggested a much more gradual sSFR-evolution than expected by most theoretical models (see, 
e.g., IWcinma nn et al. 120111 an d references therein). The most recent efforts have focused on quantifyin g the impact 



of nebular emission lines (e.g. iSchaerer fc de Barrosl I2010L Ide Barros etaLl 120121 iGonzalez et all 120121 IStark et al.l 
120121 ) on stellar mass measurements. It is presently unclear, however, whether corrections for nebular emission cause 
significant deviations fr om the nearly flat evol ution that was found prior to their implementation: line-corrected sSFR- 
values scatter about the IBouwens et al~l (|2012l ) measurement^ at z < 6 and only then become consistently larger than 
non-corrected ones. Our analytical parametrization of the sSFR-evolution does not trace the apparent increase at 
z > 6 but these high redshifts are not the main focus of this article. 
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